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Large cosmic-ray bursts as recorded with an unshielded ionization chamber have been found 


to be coincident with extensive atmospheric (Auger) showers as detected with Geiger-Miiller 
counters. The experiments show that the bursts observed in the ionization chamber represent 
regions of very high particle density in an Auger shower. The integral size-frequency distribution 


curve of these bursts is approximated by an inverse power law with an exponent 2.15+0.35. 


CARNEGIE model C ionization chamber! 

was used to detect cosmic-ray bursts of 
more than 100 particles. The chamber was com- 
pletely unshielded, and it was located in the 
greenhouse of the Botany building on the campus 
so that there was little shower-producing ma- 
terial above or to the side of the chamber. In 
order to detect extensive showers, five coincidence 
units were used, and the G-M counters for these 
were arranged on light wooden racks suspended 
over the chamber. The units, AA and BB were 
high resolving power 2-fold coincident circuits 
(1.2 and 2.0 microseconds resolving time, re- 
spectively); CCCC and DDDD were 4-fold cir- 
cuits, and E was a 5-fold circuit. A typical ar- 
rangement of the counters is shown in Fig. 1. 
Coincidences of the circuits were recorded by 
means of small lamps on the same photographic 
film as that on which the bursts were recorded. 
Thus a continuous photographic registration was 
obtained both for the bursts and for each 
coincidence circuit. Special care was taken to be 
certain that each circuit operated independently 
of all others, and the examination of the photo- 


1A, H. Compton, E. O. Wollan, and R. D. Bennett, 
Rev. Sci. Inst. 5, 415 (1934). 
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graphic record allowed an excellent check on the 
performance of the circuits. 

A total of 107 bursts were observed while the 
counters were in operation; the number of 
particles in these bursts ranged up to 3000. Of 
the 107 bursts, at least 94 were observed to be 
coincident with a simultaneous tripping of the 
counter arrangements. It should be noted that 
the counter arrangement shown (Fig. 1) was 
modified during the experiment by increasing the 
separation of counters. Consideration of the 
typical arrangement in the diagram shows that a 
shower of high particle density is necessary to trip 
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Fic. 1. Arrangement of the ionization chamber and 
G-M counter sets. The dimensions of the counters in each 
rack are as follows: Rack 1—2X19 inches; Rack 2— 
2X 15 inches; Rack 3—2 X 10 inches; Rack 4—1 X 10 inches. 
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Fic. 2. Cumulative size-frequency distribution curve for 
bursts observed in an unshielded ionization chamber. 


all of the counters. Such showers correspond to 
those investigated by Auger.? Thus the large 
bursts which occur in an ionization chamber with 
negligible shower-producing material near it 
represent a high density region of an Auger 
shower. Since the Carnegie chamber has an 
effective area of 0.1 sq. meter, these densities 
range from 1000 to 30,000 particles per sq. meter. 
It was observed that some of the Auger showers 
which tripped all of the counters did not have a 
sufficient particle density in the region of the 
ionization chamber to produce a burst in it. 
Since the ionization chamber is a unique in- 
strument for the measurement of high particle 
density in a shower, it is planned to investigate 
the high density region (the so-called “‘core’’) of 
an Auger shower by using 2 ionization chambers 
in a coincidence arrangement. Thus it should be 
possible to measure the extension of the core. 
Over a period of 2153 hours, 170 bursts were 


2P. Auger, R. Maze, and T. Grivet-Meyer, Comptes 
rendus 206, 1721 (1938). 
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observed in the chamber. The cumulative? size- 
frequency distribution curve of these bursts 


(Fig. 2) shows the frequency (J) of bursts con- 
taining more than S ionizing particles as a func- 
tion of S for values of S larger than 100. It is seen 
that on the double logarithmic plot the experi- 
mental points can be approximated by a straight 
line representing the inverse power law: 


N(>S)=C-S~’, 
where 
2.15+40.35. 


The apparent discrepancy whereby the smallest 
size bursts fall out of line with the curve will be 
discussed in a more general paper on bursts to be 
published shortly. The problem of bursts under 
35 cm of Fe and under intermediate shields of 
iron around the ionization chamber will also be 


discussed in the same paper. 
From shower theory one can estimate the 


energy that a primary electron or photon must 


have in order to initiate at the top of the atmos- 
phere a shower which is able to reach sea level 
with a high particle density in its core. Thus a 
burst of 100 particles represents a shower initi- 
ated by an electron with energy of about 10" ey, 
It is estimated that the largest observed burst 
(3000 particles) must originate from an electron 
or photon of at least 10'®-ev energy. According to 
the cascade theory of Auger showers, an electron 
initiates a shower close to the top of the atmos- 
phere; this shower then develops until it has 
penetrated to sea level, where it extends over at 
least 10,000 square meters. The frequency of the 
largest observed bursts (one in 500 hours) as 
compared with the frequency of Auger showers as 
measured with counters, shows that the very high 
density core of the showers has a relatively small 
area as compared with the total area of the 
shower. For the largest bursts, the core represents 
about 1 percent of the total area of the shower. 

The writer wishes to thank Dean A. H. 
Compton for his generous support. The problem 
was suggested by Professor Marcel Schein, and 
the writer gratefully acknowledges his continuous 


advice. 


3M. Schein and P. S. Gill, Rev. Mod. Phys. 11, 267 
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A quantitative calculation is made of the entire interaction called into play when two 
hydrogen molecules approach. The forces consist of three parts: exchange, quadrupole, and 
dispersion (van der Waals) forces. By compounding these, it is possible to account in a funda- 
mental way for the size and the shape of the molecules. Also, the interaction curves (Fig. 2) 
are in reasonable agreement with empirical curves derived from second virial coefficients. , 
Because of the complexity of the exchange force calculation, a simple state function of the 
Wang type (but with undetermined screening constant) was used in this work to represent 
the charge distribution within the molecule. Quadrupole moments are computed (Table V), 
and comments are made on the relation between the molecular problem treated in this paper, 
and the corresponding atomic problem (interaction between helium atoms). 


N the field of intermolecular forces, interest 

has chiefly centered about the long range 
attractive effects which can be calculated with 
relative ease. Except in the simplest cases of 
atomic interaction,' rather crude devices are in 
use for estimating the short range forces of 
repulsion, devices which depend in most in- 
stances on empirical knowledge of the gas kinetic 
sizes of the atoms. For molecules, no attempts 
have apparently been made to derive the quanti- 
tative aspects of the repulsive forces, despite 
the fact that knowledge of them is indispensable 
for an adequate ynderstanding of such funda- 
mental properties as molecular size, shape, 
rigidity. Nor is it possible to determine the 
position of the van der Waals minimum of the 
interaction curve without a fairly accurate 
picture of the repulsive exchange forces. 

To obtain these forces without the tedium of 
an a priori calculation, numerous workers? have 
undertaken the useful task of deriving them 
from observed data, such as second virial coeffi- 
cients, energies of crystal lattices, Joule- Thomson 
coefficients, and the like. The relative success of 
this procedure has created an attitude of satis- 
faction and apparently a waning of interest in 
the fundamental problem involved. Endeavor 


He, cf 
G, Gentile, Zeits. f. Physik 63, 708, (1930); H 
Phys. Rev. 56, 1000 (1939). 

* Among them are: J. E. Lennard-Jones, Proc. Roy. Se. 
A106, 463 (1924). For a full review of this work see R. H 


. J. C. Slater, Phys. Rev. -. 349 (1928) ; 


. Margenau, 


Fowler, Statistical Mechanics. M. Born and J. E. Mayer, 
Zeits. f. Physik 75, 1 (1932). J. O. Hirsch elder, R. B. 
Ewell, and 4 R. Roebuck, J. Chem. Phys, 6, 205 (1938). 


has been shifted toward correlation of diverse 
gas kinetic phenomena by means of a single set 
of interaction curves, without much regard to 
the basic credentials of these curves. The 
present paper is a step in the opposite direction 


‘inasmuch as it derives interaction curves from the 


elements of quantum mechanics. This may entail 
a sacrifice in numerical accuracy of the results, 
but it should fill the gap between semi-empirical 
reasoning and theoretical understanding. 

In the case of molecular hydrogen, the total 
interaction consists of three main types: First, 
the exchange forces just mentioned which owe 
their origin to the interpenetration of the elec- 
tronic clouds; they are responsible for the 
rigidity of the molecule. In the calculations they 
appear as first-order effects, arising when state 
functions of the proper symmetry are used. They 
depend very strongly on the relative orientation 
of the interacting molecules. The major part of 
this paper is devoted to their study. 

Secondly, there are the forces resulting from 
the presence of a permanent quadrupole moment 
in the Hz molecule. Their dependence on orienta- 
tion is also strong, but they vanish in the mean 
over all orientations. In comparison with the 
other types, the quadrupole forces are small 
except at distances of separation much greater 
than those of interest.in connection with gas- 
kinetic phenomena. In the calculation, they 
appear in the same formalism as do the exchange 
forces, being first-order effects and associated 
with the symmetry of the molecule. Their 
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magnitude is a sensitive function of the concen- 
tration of electronic charge about the nuclei. 

In the third place, the attractive van der 
Waals forces are to be included.* As is well 
known, they represent the instantaneous attrac- 
tions between the electronic multipoles as they 
rotate, within one molecule, in partial phase 
agreement with the multipoles in the other. 
These may be calculated by a second-order 
approximation method and with the use of state 
functions lacking the correct symmetry. In the 
H; problem, it is necessary to include in the 
treatment both dipoles and quadrupoles,’ and 
it is safe to neglect the interaction of higher 


multipoles. 


I. EXCHANGE FORCES; NOTATION AND METHOD 
OF CALCULATION 


The method used in the calculation of the 
exchange forces is based on the work of Slater ;* 
some of its features, as they relate to the four- 
electron problem, are discussed by Glasstone, 
Laidler, and Eyring.’ In order to achieve 
simplicity, these authors neglect all multiple 
exchange effects, thereby reducing the number 
of exchange integrals from 24 to 7. Unfortunately 
this curtailment invalidates all quantitative 
conclusions one might wish to draw, for the 
exchange integrals do not arrange themselves in 
descending order of magnitude as the number of 
transpositions characterizing the exchange in- 
creases. They are, in fact, nearly all of com- 
parable magnitude ; even if the non-orthogonality 
between orbitals is small, exchange integrals 
may well be large. Furthermore, in the case of 
the hydrogen molecule, the non-orthogonality 
integral has the value 0.72 and cannot be 
neglected. 

In order to carry through the work, it appears 
necessary to use the simplest possible type of 
state function for the Hz molecule. A function 
similar to that employed by Wang® was therefore 


*For a review of the theory of these forces see H. 
Margenau, Rev. Mod. Phys. 11, 1 (1939). Further calcu- 
lations dealing specifically with He were made by H. S. W. 
Massey and uckingham, Proc. Ir. Acad. 45, 31 (1938). 
See Sec. V of the esent paper. 

. . C. Slater, P - Rev. 38, 1109 (1931). 

Glasstone, K. J. Laidler, and H. Eyring, The 
Theory of Rate Processes (McGraw-Hill, 1941). 
8S. Cc Wang, Phys. Rev. 31, 579 (1928). 
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chosen. It will be seen, however, that the value 
of the nuclear charge Z, which gives the lowest 
energy for a single molecule, leads to erroneous 
consequences in the present problem. The mean- 
ing of this will be discussed in due course; we 
note at present that this fact enforces the 
preliminary use of an undetermined Z. 

Let the four protons in the two interacting 
molecules be labeled a, 5, c, d. They are first 
taken to be located at arbitrary, fixed points. 
The symbols a to d will also be used to represent 
the hydrogenic wave functions of an electron 
about the nuclei a to d. Thus, for example, 


= (Z*/mao*)! exp (—Zre/ao), (1) 


r, denoting the distance of the electron from 
nucleus b. We note that this function satisfies 


the Schrédinger equation 


h? Ze? 
— (2) 


The electron coordinates will be numbered from 
1 to 4, so that, for instance, 5(3) represents 
electron 3 centered about nucleus b. If the spin 
functions a, 8 are introduced in the usual way, 
a product function like 


v1 (3a) 


signifies a unique assignment of electrons to 
nuclei. From this product function, an anti- 
symmetric Pauli determinant may be con- 
structed; this, when normalized, will be written 


in the form: 
laa ca dB|. (3) 
The state of the four atoms under considera- 
tion is a singlet state, being composed of two 
molecules in singlet states. The only functions 
of type (3) which can cooperate in forming that 
state are those having a total spin 2,=0, and 
they are 6 in number: 


bB dp\, 
ba 
bg da|, 
da|, 
da|. 


Vi=|aa 
W2=|aa 
V;=|aa 
ba 
Ws=|a8 ba 
bp 
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’ ‘TaBie I. List of symbols used for (1) /(Pigi)Heidr and (2) S'(Peei)¢idr. Symbol for integral is given opposite the 
permutation Px (which is written in the form of cycles). Unprimed symbols refer to (1), primed symbols to (2). 


(ab) 
(ac) 


U1, Uy" 
(ad) | 
| 
| 


(abe) 


| 
| (bac) 


U2, Us’ 
U3, Us’ 
(bc) 
(bd) 
(cd) 


U4, U4 


Us, Us 


Ug, Us’ 


(abd) 


| (ab)(cd) Wi, wr’ 21, 21" 
(ac)(bd) we, we’ 


(ad) (bc) | wa, ws’ 


(abcd) 
(dcba) 


| 
| £1, 21 
(bacd) | 
| 


V1, 


Ze, 


. 


Ve, v2" 
V2, V2" 22, 


23, 33 


V3, v3" 23, 
V4, 04" 


U4, U4 


| 
| 
V3, 03" 
| 


Of these, two singlet functions can be constructed 
by known rules:’ 


=F (Vi —Vs— Vit Vs), 
Vp=}(V2—Vs— Vit Vo). 


Between V4 and Wz the Hamiltonian is to be 
diagonalized, and this process leads to the 
secular equation 


Hap—ESap 
Has—EAas 


When the matrix elements appearing in this 
equation are expanded with the use of (5), each 
of them becomes a linear combination of ele- 
ments and respectively, where the 
subscripts 7,7 refer to the set of functions (4). 
The evaluation of H;; and A;; proceeds as follows. 

From the properties of the antisymmetric 
functions and the symmetry of H it is clear that 


(5) 


=0. (6) 


(7) 
where P, is some permutation among the four 
electrons whose coordinates are contained in ¢;j, 
and p, is 1 for odd, 2 for even permutations. 
Also, ¢; is written for the simple product function 


of type (3a) corresponding to the y; of type (3), — 


and the integral includes summation over the 
spins. But an integral like {y1P,g,r can 
always be reduced to the form 


f 


with ¢; given by (3a). We see, therefore, that 
7 See, for example, reference 5. 


every H;; is a simple sum of several exchange 
integrals of this form. These fall, of course, into 
5 classes, in accordance with the properties of 
the symmetric group on four particles. It is well 
to enumerate and label them by reference to the 
permutation P, characterizing the particular 
integral in question. 

In doing so, however, the spin functions may 
be omitted from ¢;, for they either disappear in 
the summation or cause the particular integral 
to vanish. Thus, for example, we shall mean 
henceforth by 


f 
simply : 


f 
X 


so that when P;, is (abc), for example, 
f 


In Table I we list the names which will be 
given to the integrals /P,¢i/7¢idr opposite the 
P, themselves; the latter are written in the form 
of cycles. 

The equality among some of the integrals is 
due to the Hermitian character of the Hamil- 
tonian; this causes integrals belonging to a 
permutation and to its reciprocal to be equal. 
No further equalities exist unless the position 
of the nuclei possesses some symmetry, but this 
cannot be assumed as yet. 

The decomposition of H;; into exchange inte- 
grals is effected by the use of formula (7) ; many 
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TABLE II. Matrix elements H;; in terms of exchange integrals. 


| wi—22:+w; 


€— — Ue +, 


1 

2 €— U1 — Ug t+ | | | 

3 €— U3— Uy wi+we— | 
5 €—U2—Us+We 
6 


permutations are absent from the sum because 
of the orthogonality of the spin functions. In 
Table II all elements H;;(=H;;) are collected. 
When the matrix elements H4s, Has, Hes, 
which enter into the secular equation, are 
_ evaluated by means of Table II, the subsequent 
expressions result. 
+witwe 
(21+22— 223), 
— (vi +02+03+04) (8) 
+ws— (22-+23— 221), 
+wit+w2 
+w3+2(21+23— 222). 
Attention must also be given to Aga, Aug, and 
Azs. But these are simply related to the corre- 
sponding H elements. If we define quantities 
u;’, w;’, 2;/ to be identical with the un- 
primed quantities given in Table I except that the 
operator H in the exchange integral is omitted, 
then the A elements result from the 7 elements 
on priming all quantities on the right of Eqs. (8). 


Il. EXCHANGE INTEGRALS 
The exchange integrals ¢, “;-- +23 listed in Table I may be decomposed into elementary exchange 


To proceed further, it is expedient to adopt 
some conventions regarding the detailed: treat- 
ment of the Hamiltonian operator in the four- 
electron problem. 

In writing the Hamiltonian we use the follow- 
ing notation: Let ra; be the scalar distance 
between the ith electron and nucleus a, and put 
(9) 


Then 
2 4 


H=-—)>V? 


2m 


‘ 
1 


i>j 


provided Ey stands for the repulsive Coulomb 
energy between the four nuclei. With the use of 
Eqs. (1) and (2) we now obtain 


Z*e? 
Ha(1)b(2)c(3)d(4) = 
0 


—(a2tastay 
vit 61+ 6245s) 


+E put Ew (10) 


i>j 


integrals, which will now be considered. While in the former the integrand is the entire Hamiltonian 


(10), the latter contain only one of its terms. To label them both succinctly and naturally, it is 


expedient to use the abbreviations (9). Thus we shall define 


(ayd) = 
el 


2 
12 


= f a(1)b(1)drs, 
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and others in a similar manner. As to symmetry, it will be noted that 
(ayd) = (dya), 
(aab) = (a8b), 
(abped) = (chpad) = (adpcb) = (cdpab) = (bapdc). 


That is to say, the order of the Roman letters on one side of the Greek letter may not be changed 
if that on the other remains unaltered, but a Roman letter may be shifted freely from one side of 
the Greek letter to the other, provided its order is unchanged. 

The elementary exchange integrals fall into four categories, in accordance with the number of 
nuclei whose positions affect their value. All one-center integrals, like (aaa), (b8b) etc., are of course 
equal and independent of the position of the nuclei. Among the two-center integrals one may dis- 
tinguish four different types, exemplified by the following specimens: 


(aBa), (a8b), (abpab), (abpba). 


These all occur in the problem of the hydrogen molecule and are well known. 
Three-center integrals fall into three classes, characterized by 


(afc), (acpbc), (acpcb). 


For special positions of the nuclei, they have been calculated by Gordadse,* Hirschfelder, Eyring, 
and Rosen,® and by Coulson.” Since they are of rather complicated structure and, in their exact 
forms, quite unsuggestive, it is important to have a quick way of estimating their magnitude. 
Consider, for example, (a8c). The product function a(1)c(1) is largest in the region halfway between 
the nuclei a and c; therefore, if 6 is far from a@ and c, 


(afc) ~ Aace®/ Rs, acy 
where Rp, ac is the distance from 6 to the midpoint between a and c. In a similar way, 
(acpbc) ~ (ab) ~ Aase®/ Re, ab, 
(acpcb) ~ Rac, 


These approximations are quite accurate when the distances R involved in the formulas are several 
times as large as the arguments of the A functions. The integrals do not, however, become infinite 
as R-0; they take simple limiting forms which are either two-center integrals or else very manageable 
expressions easily obtainable from references 8 and 9. 

For more accurate work, approximation (11) is not adequate. On the other hand, the general 
forms of (acpbc) and (acpcb) which are needed in this work are not known. Their exact calculation 
would seem to entail rather formidable labor. Fortunately, however, it is possible to avoid most 
of it by a simple reduction process now to be described. 

If the function c* is contracted more and more about its nucleus, until finally it becomes a 6 function 
located at the nucleus, the integral (acpbc) turns into (ayb). Alternately, if the function ab is approxi- 
mated by a 6 function located midway between nuclei a and 6, the same integral reduces to the 
two-center integral" (caBc)-Aas. When one distance is larger than the others—as is true in the 
interaction of two molecules—these reductions can always be made and lead to nearly equal numerical 
results. This we regard as a test of the validity of the reductions, and we take 


(acpbc)=3[(ayb) J. 


(acpcb) =} [Aae(caryb) +Aa(ayBc) J. 


*G. S. Gordadse, Zeits. f. Physik 96, 542 (1935). 
* J. Hirschfelder, H. Eyring, and N. Rosen, J. Chem. Phys. 4, 121 (1936). 
”C. A. Coulson, Proc. Camb. Phil. Soc. 33, 104 (1937). 


1 @B here stands for the reciprocal of the electron distance from the midpoint between a and b. 


(11) 


Similar reasoning shows that 
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The remaining integrals are known.’ To give an example of the extent to which this procedure is 
reliable we consider the case in which the nuclei a, 6, c are collinear, c is between a and b and the 
distance c—b is three times the distance c—b, the latter being the internuclear distance in He. The 
two terms in the above expansion of (acpcb) then have the values 0.1448 and 0.1417 atomic unit, 
respectively. For larger distances ¢—b, the reduction works still better. 

All four-center exchange integrals are of the form (abpcd). Their exact calculation would be 
extremely tedious if it can be carried out at all. But here again the scheme just described will be 
used. The reduction can be effected in two ways, and if the results do not differ beyond the tolerated 
limit of error, their mean is taken as the value of the four-center integral. Thus 


(abped) ~ }[Aac(baryd) + Ava(aBic) 
When a and c belong to the same molecule, the two terms on the right become very nearly equal, 
and each approaches A,r*(e?/R) in the limit in which R, the distance between molecular centers, 


is large. 

We do not wish to lengthen the present paper by an inclusion of the detailed formulas for the 
residual integrals, which may, if they are found useful in other molecular problems, form the sub- 
stance of a later communication. As far as our present purposes are concerned, the reader will be 
relieved to know that an even simpler scheme than the alternate reductions here studied, a scheme 
which will be described in a later section, leads to significant results for the interaction energy of 


molecular hydrogen. 
When the exchange integrals ¢, u;---zs3 are expressed in terms of their elements, certain combi- 


nations often occur. These will first be singled out. We define 
A = (aaa) = Ze?/ao, 
D=—}3Z%e?/do, 
Ba=e?/Ra— (aa), 
Ca» = Rav + (abpab) — 2(aBa), 
Tar = (a8b)/Aaw, 
(abpba) /Ags— 2T os. 
Here C,» will be recognized as the ‘‘Coulomb energy,” Xa» as the ‘‘exchange energy’’ between two 


hydrogenic atoms at a and b. All these expressions are functions of only R.», the distance between 
a and b, and most of them vanish exponentially at large distances. In addition we need the following 


more complicated combinations : 
X= 2[(acpbe) + (adpbd) — (ayb) — (a6b) Aus, 


e 2 


+ (abpbc) /AesAse+ (acpba)/AarAac+ (bcpca)/ 
Aw Arve Mac 
2 (acpbd) 


(4 e 
te te 2[(ayb) + (a6b) ]/ Aas — 2[ (cad) + (cBd) Aca, 


(bcpab) (bdpac)  (bapad) (bepcd) (capbd) (dapcd) 
—[(ayb) + J/ — [ (bac) + (b5c) ]/ Ave — [ (cad) + (cBd) Bd) + (ayd) }/ Aaa. 


Watea= En + 
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By permutation of subscripts, other functions can be constructed from this list. Not all of these 

are different; for instance it will be seen on inspection that there are four different U functions 
(Wore Uaca, Urea), three different V functions ( Vo», ca, Vac,ra, and three different W 
functions (Wavea, Weaca, Waera). Only X’ and U’ have the property of vanishing when a reduction 
of the type (11) is performed, which may be taken as an indication that they are generally small. 
They disappear when the nuclei in each molecule coalesce. 

The function V.»,-¢ has a rather interesting significance. When it is ‘‘reduced,”’ it represents the 
electrostatic interaction between two linear quadrupoles, one consisting of protons at a and 6 with 
two electrons at the midpoint between them, the other of protons at ¢ and d and two electrons at 
their midpoint. We shall show later that this is indeed its true significance, and that the value of 
this integral is well approximated by this quadrupole interaction. 

With the use of these abbreviations, the exchange integrals take on symmetrical—even if some- 
what complicated—forms, viz. : 


€=4D+4(Z—1)A + Cat Cact Coat Coat 


= [4D+2(Z—1)(A + Tac) +X 
us= Asal 4D+2(Z—1)(A +T aa) +X 
a= 4D+2(Z—1)(A Bas + Bact Brat+Beat+ 
4D +2(Z—1)(A + Tra) +X Boa +Bre+ Beat 
4D+2(Z—1)(A ca) +X cat Boc+ Boa cal, 


01 = 4D+Bea+ Bra t+Bea—(2—Z)(Tas+ Tre) +(Z—1)A + Uare+ Ure], 
02 = 4D + Bact Bre+ Bea—(2—Z)(Tav+ Toa) +(Z—1)A + Vesa], 
4D +Bas+Bre+ Boa—(2—Z)(Tact+ Tea + Toa) +(Z—1)A + Uscat+ Urea 
4D+ Bas Tea t+ Tra) +(Z—1)A + Usea+ Urea], 


(14) 


2(Z— Tea) + Vad, ca], 
1)(Tae+ Toa) + Vac, ba], 
4D +X eat Xve+2(Z— 1)(Taa+ Tre) + Vea, re], 


31 = 4D — (2 —Z)(Tast+ Tre+ Tcat+ + Wares}, 
4D — (2 —Z) (Tact Tar + Tra) + Woaca 
23 = AacAsaArcAaal 4D — (2—Z) (Tact Treat Tre+ Toa) + Wacra ]. 


Thus far our development has been quite general. We can now afford to make specific assumptions: 
Protons a and b belong to one molecule, protons c and d to the other; the distances a—b and c-d 
are equal. Functions with subscripts ab and cd may therefore be written without subscripts, the 
understanding being that they are to be evaluated for the internuclear distance of the Hz molecule. 
It is also possible to distinguish different orders of magnitude among the various constituents of 
Eqs. (14). All quantities without subscripts (or with subscripts ab, cd) are large, all others small. 
Whether a term can be totally neglected can only be decided by inspection of the A functions which 
multiply it.’ This overlap integral has the well-known simple form 


Aas=A(s) = (1+s+4s*)e“, 
Ss being (Z/ao)Ras. The ratio Ag-/Aas is therefore always fairly small, and it is safe to neglect its 
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fourth power, though not the second. We are thus enabled to simplify the list (14) somewhat, as 
follows: 


Coat Coa, 


Boa], 
y= xe], 

oa], 

y= 


01 = + + Usne— (2—Z) (Tact Tre) 
AAsaAaal 4D +B —(2—Z)T+(Z—1)A + Ussa—(2—Z) (Toa + Toa) 
03= 4D +B —(2—Z)T+(Z—1)A + + Unca—(2—Z) (Tact Tad) 
04= Ady Aral 4D +B—(2—Z)T+(Z—1)A + Usea t+ Tea) 


(14b) 


W2=w;=0, 


= A*ApAaal 4D — 2(2—Z)T — (2—Z) (Tre +T aa) + Warea ], 
—2(2—Z)T —(2—Z) (Tact Tra) + Weaca 


23=0. 


In “; the term X’ has been dropped because of its smallness, and in u2 to us, the C functions are 
neglected. Numerical comparison with terms retained will justify this curtailment. 

The primed functions, defined in Table I, are easily obtained from this list; for they are simply 
the product of A functions which appears in the equations for their unprimed mates. Thus, for 
example, e’=1, u;’=A?, v»4=AA,-Aza etc. We are now ready to compute the two roots E of Eq. (6). 


Ill. GENERAL EXPRESSION FOR EXCHANGE ENERGY 


Of the two solutions of the determinantal Eq. (6), the one of interest is that which will represent 
the energy of two H:2 molecules at infinite internuclear distance. It is convenient, therefore, to 
solve the equation first of all ‘‘in zeroth order,” i.e., with the neglect of all terms which vanish at 


infinite separation. 
As a preliminary, let us recall the results of the Heitler-London treatment of the single molecule 


and translate them into the present notation. Here we have 
| aa bB|, |ap ba|, 
and Wa =yYityo. The function y, represents the-triplet, y_ the singlet or normal state of the molecule. 
In view of Eqs. (2) and (12) 
Hab =[2D+(Z—1)(a1+ 82) jab 


f 


so that 


whereas 
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For the ratio of these, which represents the energy of the molecule, we have 
E(A2) =2D+ (1 (15) 


where the lower signs refer to the stable state. The second part of this expression represents the 
molecular energy since 2D is the energy of two hydrogen atoms of nuclear charge Ze (cf. Eq. (12)). 
Wang has found the value of Z which minimizes (15) to be 1.166. 

We now turn to Eq. (6), in which the H elements must be expressed in terms of the functions 
(12) and (13) via Eqs. (8) and (14). Before doing this it is well to note the effect of subtracting 
Eda from Haa. From what has been said about the A elements it is apparent that this subtraction 
merely amounts to the replacement, in H4a, of every term 4D by 4D—E. The same is true regarding 
Has—EAas and Hpg—EAgs. We therefore define new matrix elements, Kaa, Kas, Kes, which 
differ from the corresponding components of H merely in the substitution of 


—E’=4D-E 
for 4D. According to its definition, E’ is the total molecular energy of the two interacting partners, 
including the binding energy of each molecule. Equation (6) now reads 

Kaa Kas 

Kas Kap 


In zeroth approximation, only those parts of (14b) need be retained which bear no subscripts. 
This leaves us only with parts of ¢€, “1, %s, and w;. When these are used in (8), we find 


—E'+2C+4(Z—1)A X+C+2(Z—1)(A +7) ] —E'4+2X +4(Z-1)T], 
Kaa=}Kua. 
On insertion of these expressions into (16) there results ; 

=0, 


=0. (16) 


and this allows the two roots to be determined by putting 
=0 and ius. 


When these are solved for E’, the first leads to 2E(H2) as given by Eq. (15) with the choice of the 
positive signs, the second to the same equation, but with negative signs. Our interest is therefore 
confined to that root of Eq. (16) which is given in zeroth approximation by the equation 


Kia=0. (18) 


We now show that the interaction energy which is being sought is similarly a solution of 
Kaa=0 (19) 


so that we need not calculate the elements Kzz and K4, at all for the purposes of the present problem. 


Let us put 
Kaa=—LE'+Lo+lE’ +l 


where —LE’+Lp> represents ls (the coefficients L and Ly are easily identifiable from (17)), the re- 
mainder being small. We may then treat L and Ly as “‘large,” J and J) as small coefficients. Similarly 


Kas= — ME'+Mo+mE'+m, 
and we know from (17) that M=4L, My=4Lo. Finally, 
Kzgs= — NE'+No+nE'+no. 
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In solving (16), it is now convenient to write 


E,’ being the solution of (18), i.e. the energy of the two Hz molecules at infinite separation. The 


determinant (16) then reads 
— LEo' +Lo—Ln+lEy’ +lo 
— ME,’ +Mo—Mn+mEy' +m 


— ME,’ + Mo—Mn+mE,'+mo 0 


In expanding it, the squares of small terms may be omitted, and we know that only the first two 
summands of each element are large. Also, Eo’=Lo/L. Insertion of this causes all large terms in 
the elements of the first row to disappear because L/L)»>=M/Mpo. Expansion then gives 


+h), (20) 
and this is simply the condition (19). We turn, therefore, to the evaluation of K44. When the terms 


are suitably grouped, there results 


Kaa=— (144%) +[$E0—B—(Z—1)A Ava) 
+ Coat Cro+ Cra + 2A*(Bac+ Boe + Bre+ Bra) — ae] 
Toa] — Tre] 
— Xsa+Xre+2(Z—1) Tra + + 
Uare+ Unse— (2 -Z)(Tac+T ve) ] 
—AAsaAval Uasa+ Ussa— (2 —Z)(Taa+ Tra) ]— AAacAaal Vaca + Urea —(2-Z)(Tac+ Toa) 
— AdscAsal Usea+ Usea— (2 Toa) Aachsa) [Eo+2(2-Z)T] 
—A*Ap-Aaal, Wabea— (2 —Z) (Te+ ]—A*Aachval Woaca — (2 —Z)(Tac+ Toa) J+A‘ Vas,ca- (21) 


n is the exchange energy we are seeking; it is found by equating K4a to zero. 


IV. NUMERICAL RESULTS 


When preliminary computations based on the 
use of Eq. (21) were made, a rather curious fact 
emerged. At first Z was taken to be 1, then 1.166 
(Wang’s value), and in both cases did » become 
negative at all significant distances of interaction. 
In this work, alternate reductions were made as 
described earlier, and careful attention was paid 
to the approximations involved; the conclusion 
was reached that the negative result was not 
occasioned by inaccuracies in the numerical 
evaluation of the integrals. To test the matter 
further, the integrals constituting » were simpli- 
fied by allowing the nuclei of each molecule to 
coalesce. The result was then identical formally 
with the interaction energy of two helium atoms. 


In that problem, too, 7 is negative for Z=1 and. 


1.166. It is possible to show, moreover, that the 
helium interaction, calculated by means of 
hydrogenic functions, is repulsive only for values 
of Z above the critical value Z=11/8. Since 


this appears to be a somewhat interesting feature 
of the 4-electron interaction, we shall prove it in 
a later section. For He, of course, the matter is 
unimportant since the value of Z which mini- 
mizes the atomic energy is 27/16, well above the 
critical value. One would expect, then, that in 
our problem Z also possesses such a critical 
value, presumably in the neighborhood of 11/8. 

In a sense, this is a rather sad commentary on 
the adequacy of the Wang function for the 
purpose at hand,” which might induce one to 
look for better wave functions. The complexity 
of the present calculation, however, definitely 
counsels against that undertaking. The situation 
appears even more embarrassing when it is noted 
that Wang’s value of Z describes the polariza- 


2 A more favorable, literal interpretation of this result 
would entail the consequence that hydrogen molecules, 
while indeed attracting each other dynamically at nearly 
all distances, are kept apart for reasons of entropy. This, 
however, we do not believe to be the case. 
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bility of Hz surprisingly well," gives a good 
account of its magnetic susceptibility and a 
reasonable one of the dispersion forces." 

But on further reflection, the puzzling aspects 
fade away. The effective nuclear charge has no 
definite physical significance—except insofar as 
it measures roughly the concentration of the 
electron cloud—and may well take on different 
values in different physical problems. In partic- 
ular, there is nothing unique about the Z which 
minimizes the molecular energy. It happens that 
in the calculations of the polarizability, the 
susceptibility, and the dispersion forces the mean 
square of the charge distribution is of importance, 
and this would account for a common Z. The 
present problem, however, has altogether differ- 
ent features. 

The change in Z from 1 to 1.166 which occurs 
as we pass from the free hydrogen atom to the 
molecule is a measure of the increased electron 
concentration about each proton occurring as a 
result of the repulsion by the other electron. 
The term “screening,” though quite adequate 
in atomic problems, loses its significance here. 
Now, when two molecules interact, there is not 
only one electron which tends to drive a given 
electron back to its nucleus, but three. This 
crude picture may be made to give a qualitative 
indication of the value of Z to be chosen in the 
present problem: If one electron causes an 
increase in Z from 1 to 1.166, three might 
increase it from 1 to about 1.5. A similar result 
could be argued in this way: as the two nuclei 
of the Hz molecule coalesce, Z changes from 1 
to 2; as the four nuclei here under consideration 
are made to coalesce, Z changes from 1 to 4 
(except for screening, which is neglected in both 
instances). Wang’s value of Z represents a 
stage in the process of coalescence 1/6 on its 
way toward completion ; the corresponding stage 
in the four-electron problem is given by Z=1.5. 
This result is likely to be too large, of course, 
because the distance between nuclei in the 
different molecules is greater than in H2. 


¥J. O. Hirschfelder, J. Chem. Phys. 3, 555 (1935). 


We have computed the polarizability of Hz as a function 
of Z and have found that Z=1.167 gives about the best 
fit with observation. 

“H.S. W. Massey and R. Buckingham, Proc. Ir. Acad. 
45, 31 (1938). 
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While we do not wish to ascribe much quanti- 
tative significance to the present argument, we 
are forced, in the absence of more conclusive 
evidence, to use it. The exchange forces do not 
depend as sensitively on Z as, for example, the 
second-order forces. Our results will not be 
changed decisively by a change in Z from 1.4 to 
1.5. In the following we select one approximately 
“right” value, Z= 1.428, and one “‘wrong”’ value, 
Z=1.785 for comparison. Their exact choice was 
dictated by numerical convenience; the first 
makes the parameter Zd/ao, wherein d is the 
internuclear distance of He, equal to 2, the 


~ second to 2.5. 


The preliminary calculations further indicate 
that a certain simplification is permitted. Among 
the functions occurring in 7, there are many 
which vary slowly as the nuclei of each molecule 
are made to coalesce. In these we may proceed 
to the limit in which Ra=R,.z=0. Then, if we 
denote the distance between molecular centers 
by R, 


all X’—0, 
U’—0, 
(R) +B(R)+M(R) 
Waea—>X (R) +2B(R) 
+C(R)+2M(R)—2T(R) 


The other three U functions are formed by 
fairly obvious permutations of subscripts. The 
new quantity M which appears here is given by 


Ma.= 2(aapac)/ Ace. 


M(R) may be approximated by either 27(R) or, 
better, by 7(R)+(2e?/R)—B(R/2). We finally 
introduce the abbreviations E(abc) for the 
electrostatic repulsive energy between protons 
a, b, and c (and similarly E(abcd) for the former 
Ey), 

x=Z-1 
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TABLE III. Values of numerical constants used. All 
energies are stated in units Ze*/ao. 


1.428 1.785 
Zd/ao 2.0 2.5 

A 0.5865 0.4583 
A 1 1 

B 0.0275 0.0094 
Cc — .0191 — .0128 
xX — .349 — .348 
T 0.692 0.626 
Ey’ 1.370 2.794 


and obtain from (21) 

n= Boa) 
+ Cact Coo+ Coa 
+ Aca Ava) 
Aachaa+ AscAsa) 
— B(R) — 2xT(R) + (2e?/R) ] 
— A[Aa-As-E (abc) + (abd) 
+ (acd) + Av-AsaE (bed) 
t+ 


+2(1—«)T—X (R)—2B(R)—C(R) 
2 


—2«T(R) — (abcd) (22) 


For reference, we recall : 
C+A2X +2«(A +A?T) 
1+A? 


The numerical results now to be reported were 
based on this approximation. Equation (22) is 
slightly less sensitive to orientation of the mole- 
cules than (21), but the economy of labor which 
it effects recommends its use despite this fault. 
The values of constants are collected in Table 
III. From the last entry the binding energy of 
the molecule for the assumed value of Z may be 
obtained, since, as may easily be verified, it is 
given by E,’Z/2—(Z*—1). This is —1.65 ev 
for the smaller, +8.33 ev for the larger Z. In 
the latter case, therefore, the hydrogen molecule 
would not be stable at all. 


E,’=2 


Three relative orientations of the molecules 
have been treated. In 


Case (a): both molecular axes are parallel to R, 
the line joining the centers; 

Case (b): one axis is parallel, the other is perpen- 
dicular to R; 

Case (c): both axes are perpendicular to R. 


The results for n, exclusive of the term A‘V 4s, .4, 
are given in Table IV. It is perhaps worth 
mentioning that in the region here of interest » 
can be approximated fairly well by an expression 


of the form 
n=Ae~*s 


in which 6 differs but little in all 6 cases, having 
the value 1.55+0.10, while A varies markedly 
from instance to instance. The quantity § 
represents the distance R measured in units ao/Z. 


V. QUADRUPOLE FORCES 


We have already interpreted the last term of 
Eq. (22), [A?/(1+A?) a, as the quadrupole 
energy on the evidence that V “reduces” to 
the electrostatic interaction between the two 
structures of Fig. 1 when the approximation 
(11) is made. In this approximation, then, V is 
given by the formula'® appropriate for linear 
quadrupoles, 


3 2 


in which Q, the quadrupole moment, is 2e?(d/2)?, 
and f is the function characteristic of this type 
of interaction. Explicitly, if 6; and 62 are the 
angles which the molecular axes make with R, 
the line between centers, and ¢ is the difference 
between azimuths in a plane at right angles to R, 


then 


f= 1—5 cos? 6,-—5 cos? 62.— 15 cos? 6, cos? 02 
+2(4 cos 6: cos 62+sin 6; sin 62 cos (24) 


Our calculation shows that this interaction 
between point charges must be multiplied by the 
“diffuseness factor” [A?/(1+A*)? in order to 
represent the true interaction between the charge 


clouds. 
The present approximation is in fact suff- 


15 See review article, reference 3. 
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ciently good for our purposes, but it is interesting 
to point out that an accurate calculation merely 
causes the term Q in (23) to be multiplied by 4/5. 
A calculation of this kind has already been made 
by Massey and Buckingham,“ whose work was 
limited to one special value of Z. Since their 
result is written in a rather cumbersome form 
which veils its simplicity, and our derivation is 
short, we shall give it here. ; 
Equation (23) is valid in general, provided we 
replace Q by its quantum mechanical equivalent : 


= ed<r? 3257) w. 


Here r; and 2; refer to the positions of the various 
charges e;, and the average is a quantum one. 
The sum extends over the 4 particles composing 
the molecule. We take the Z axis along the 
molecule. The protons will then contribute to 
this sum the amount —4e(d/2)*. For the elec- 
trons, (72) = (2x (a? w; hence their contri- 
bution to the sum is —4e((x*)—(z*),). The 
subscript may here be dropped because the 
electrons are indistinguishable. In view of these 
simple facts, 


d? 

The mean values appearing here are not difficult 
to calculate when hydrogenic wave functions 


are used—they also occur in the theory of 
polarizabilities.’ In general, 


TABLE IV. Exchange energy as function of 
intermolecular (R). 


s R(inA)| Pos.¢ 
Z=1.428 
6 2.22 0.358 0.249 0.150 
7 2.59 0.102 0.0602 0.0344 
& 2.96 0.0231 0.0135 0.00744 
9 3.33 0.00504 0.00287 0.00153 
10 3.70 0.00104 0.000582 0.000302 
11 4.07 0.000200 0.000114 0.0000579 
Z=1.785 
6 1.78 1.390 0.754 0.325 
7 2.07 0.426 0.192 0.0764 
& 2.37 0.100 0.0449 0.0173 
9 2.66 0.0247 0.00986 0.00354 
10 2.96 0.00545 0.00205 0.000711 
11 3.25 0.00113 0.000451 0.000137 


Fic. 1. Molecular structure and arrangement for 
calculation of electrostatic energy. 


f Pabdr | 


Here (2¢*)w= (x0?) are averages taken over the 
electron in the H atom, é and ¢ are, respectively, 
the x and z coordinates measured from the 
— of the molecule. For the functions here 


do 


1 


Remembering the form of A (cf. appendix), we 
find 


d? 


and finally 


d 
(25) 


The quadrupole moment depends on Z only 
through the overlap integral A, and in this 
simple way; the factor A?/(1+A*) appears auto- 
matically in the calculation, justifying our 
former remarks. Collecting the results of the 
present section, we write for the quadrupole 


energy 


3/ \*ed* 
i+A? 62, ¢). (26) 


The function f (cf. Eq. (24)) causes Eg to vanish 
when a mean value is formed over all orientations. 
The quadrupole moment is listed as a function 
of Z in Table V, where — Q’ is stated in units ea,’. 
The values of Eg, given in Table VI for cases 
(a), (b), and (c), are seen to be small in com- 
parison with the exchange energy in Table IV. 
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VI. LONG RANGE VAN DER WAALS FORCES 


The dipole-dipole part of the long range forces 
has been calculated for the present problem by 
Massey and Buckingham using Z= 1.166. As to 
the dependence of these forces on orientation, 
their result indicates a variation of about 50 
percent as we pass from case (c) to case (a). 
The corresponding variation in » is measured by 
a factor 3. When compounding the various 
components of the total interaction, we shall 
therefore permit ourselves to neglect the de- 
pendence on orientation of the long range 
forces.'® The numerical result for the mean over 
all orientations given by the aforementioned 


authors is 


(27) 


= 
ao R 


they also point out that it agrees reasonably 
well with empirical facts as far as they are 
available. 

In using this result in the present calculation 
we should be guilty of employing different values 
of Z in different parts of the problem. As has 
been noted, there is nothing basically contra- 
dictory in such a procedure because a variation 


TABLE V. Quadrupole moment of Hz (in units eao*) as a 
function of Z. 


1.785 


1.166 428 
0.137 


Z 1 1 
0.252 0.200 


0.284 


TaBLe VI. Quadrupole energy (Eg) as function of 
intermolecular separation (R). 


(in volts) 
Pos. @ Fo fin b Pos. ¢ 


N 


0.000361 
0.000167 
0.0000857 
0.0000475 
0.0000281 
0.0000174 


0.000504 
0.000236 
0.000120 
0.0000676 
0.0000396 
0.0000248 


— .000481 
— .000223 
— .000114 
— .0000634 
— .0000374 
— .0000232 


— .000672 
— .000316 
— .000160 
— .0000901 
— .0000528 
— .0000331 


0.000962 
0.000445 
0.000228 
0.000127 
0.0000748 
0.0000465 


0.00134 
0.000632 
0.000321 
0.000180 
0.000106 
0.0000662 


—— 
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in the best choice of Z may well be expected to 
occur as we pass from one calculation to another, 
Nevertheless we do not wish to leave this point 
without further investigation, which will also 
throw some light on the accuracy of Eq. (27). 

In the first place, it is easy to obtain Ey.4.~ for 
other values of Z by use of the variation method 
of reference 14. Taking Z= 1.428, we find 


which is almost certainly too low. Roughly, 
the dipole-dipole energy is proportional to Z-%, 
Perhaps the most accurate semi-empirical 
method? of computing the term here in question 
is to employ dispersion f values which, for H,, 
are known with considerable accuracy. For this 
purpose, London’s formula may be used : 


Sif; 


where the f’s are oscillator strengths and the 
AE’s the corresponding excitation frequencies. 
For H2, Wolf and Herzfeld!’ find that the 
dispersion curve can be represented with re- 
markable accuracy by a two-term formula in 
which 
fi:=0.69, E,=2.71X10-" erg=0.630¢?/ap, 
f2=0.84, E,=2.12K10-" erg=0.492¢e?/ap. 


These data yield 


Beaw=—10.9~ (28) 
ao 

in place of (27). This, incidentally, confirms the 
author’s former conjecture that the variational 
Hassé method, employed by Massey and Buck- 
ingham, usually gives too large an answer for 
the van der Waals force. (It does not, of course, 
contradict the variation principle!) 

But the difference between (28) and the 
empirically better substantiated result (27) is to 
be accounted for by the inclusion of the dipole- 
quadrupole interaction. A method for estimating 
the latter has been given.'* It involves the use 
of a one-term dispersion formula. To adapt it to 


16 This, in a sense, counteracts the simplification made 
in the exchange integrals leading to Eq. (22). 


17K, L. Wolf and K. F. Herzfeld, Handbuch der Physik 


(1928), Vol. 20. 
18H. Margenau, J. Chem. Phys. 6, 896 (1938). 
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FORCES BETWEEN HYDROGEN MOLECULES 


Fic. 2. Interaction energies when the approximately 
“right” value of Z (Z =1.428) is used. 


the problem at hand, we shall take mean values 
of the f’s and AE’s in Herzfeld and Wolf's 
expression obtaining the contribution 


—31X10-7*/R® erg. 


When this is added to (28), and both terms are 
expressed in ev, we have 


38 19.4 


volts, (29) 


provided R is measured in A. At R=3A, the 
d.-q. term contributes about 3, at R=2.5A 
about } as much as the d.-d. term. If, U.erefore, 
we wished to approximate (29) by means of a 
single term proportional to R~ as is often done, 
we should in this range of R (where the van der 
Waals minimum occurs) increase the d.-d. 
contribution by about 40 percent, thus obtaining 


9.0 15.4 e? 


This is not far from the result (27) which was 
obtained by Massey and Buckingham for the 
d.-d. interaction alone. 


VII. SUMMARY OF RESULTS 


In Figs. 2 and 3 we have compounded the 
partial interaction energies represented by Eqs. 
(22), (26), and (29) or (29a), obtaining a graph 
which expresses the total repulsive-attractive 
behavior of the molecules. Figure 2, drawn for 
what we have determined to be the approxi- 
mately ‘‘right’’ value of Z, should be about 


Fic. 3. Curves similar to those in Fig. 2 but with the 
“wrong” value of Z (Z =1.785). 


correct. Figure 3, on the other hand, based on 
the ‘‘wrong”’ value of Z, is presented to show 
the effect of this choice. In view of the uncer- 
tainty in this parameter, it is difficult to estimate 
the accuracy of our exchange force calculation, 
but we believe it to be not greatly inferior to the 
accuracy with which the long range van der 
Waals forces can be obtained. 

The shape of the molecule is reflected in the 
relative positions of the three minima, the one 
implying end-on collisions being much farther 
out than that for broadside collisions. Interesting 
also is the fact that the minima are of entirely 
different depths. In actual impact problems, an 
average of the three behaviors here depicted will 
have to be considered. 

Comparison with experimental facts can at 
present only be made in an indirect way. 
Lennard-Jones and his collaborators’® have ana- 
lyzed available measurements of the second 
virial coefficient of H2. The procedure is based 
on the assumption that the interaction energy 
can be represented sufficiently well by a single 
function of spherical symmetry and of the form 
(v/R*)—(u/R*). It is found that the data can 
be accommodated by several values of the 
exponent m. The curves obtained have minima 
in the region from 3.3 to 3.5A, and depths at 
minimum ranging from 2.7 to 2 millivolts. If 
curve b of Fig. 2 is taken to be illustrative of 
the ‘‘mean’’ behavior, the agreement with the 
results of Lennard-Jones is rather satisfactory. 


19See R. H. Fowler and E. A. Guggenheim, Statistical 
Thermodynamics (Cambridge University Press, 1939). 
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VIII. REMARKS ON THE EXCHANGE INTERACTION 
OF TWO HELIUM ATOMS 


The exchange forces between helium atoms 
have been the subject of several investigations ;! 
in fact the Wang function used in this work has 
been applied to the helium problem by Gentile. 
We are here interested in only one feature of 
this problem which appears to have received no 
attention, namely the dependence of the forces 
on the assumed value of Z. 

Instead of using published results?® (which 
seem to be somewhat disfigured by typographical 
errors) it will be found simpler for our purposes 
to start anew and use our compact notation. 
The ground state of a single helium atom is 
represented by 


Yo=|aa asl. 
This leads, by simple steps, to the energy 
Eo= Ho =2D—2(2—Z)A+(aapaa). 
If now we note that 
D=—-}32Z*, A=Z, (aapaa)=§Z, (31) 


all in units e’ao, then (30) will become a minimum 
for Z= 27/16, a fact which is well known. 

In the interaction problem, the ground state 
of the two atoms is described by 


(30) 


Y=l|aa aB ca cpl. 
Here the total energy becomes 
E=H/« 
with 
Al = (aacc | H | aacc) — 2(caac | H | aacc) 
+ (ccaa | H | aacc) 
e=(1—A?)?, 


A= f acdr. 
Thus we find 


—2A°[B+C+X+2(2—Z)A —2(aapaa) 
—2(3—Z)T+4(aapac)/A] 
+A‘[4(2—Z)A —2(aapaa) +4X 

—4(2—Z)T+2(acpca)/A?] 


G. Gentile,"Zeits.’f. Physik 63, 795 (1930). 


(32) 


HENRY MARGENAU 


where all functions without subscripts now refer 
to the interatomic distance R, and Eo is given 
by (30). 

With the neglect of terms in A‘, this expression 
could also have been obtained—though with 
more labor—by letting the distance s in Eq. (21) 
go to zero. Now it is permissible to neglect A‘ 
and hence the last row of (32). Since the function 
C is always much smaller and decreases faster 
than A?, the contribution to the exchange energy 
comes almost entirely from the second row of 
(32). The only terms inside the bracket multi- 
plying A? which do not vanish as R becomes 


large are 
2(2—Z)A —2(aapaa). 


If, therefore, we wish E—2E, to be positive, 
we must require 


(aapaa) >(2—Z)A. 


In view of Eqs. (31) this implies Z>11/8. 
The mathematical situation encountered in the 
H, problem when Z was chosen to be 1.166 is 
thus explained. Whether this phenomenon has 
any physical interest, as for example in the 
interaction between two negative ions of atomic 
hydrogen where Z is certainly smaller than the 
limiting value, is perhaps difficult to say. 
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X (IN UNITS 


Fic. 4. Graph for the determination of X in terms of S. 


The author wishes to acknowledge gratefully 
the computational help given him by two un- 
dergraduate students, C. E. Hummel and R.A. 
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APPENDIX 


A list of the more important functions en- 
countered in this paper is here appended. 


S=ZR/ao, 
= 


ao 8 4 6 


/ (14545). 
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X is a more complicated function involving 
Sugiura’s integral. In terms of the functions 
used and defined in Pauling and Wilson,?! 


_ 1 2K(S). K’(S) 


A(S) A%XS) 


For calculations, the graph of X given in 
Fig. 4 will be found more convenient. 


“LL. Pauling and E. B. Wilson, fr., Introduction to 


Quantum Mechanics (McGraw-Hill Book Company, New 
York). See p. 342 et seq. 
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Radioelement Sc** 
Quantum energies (Mev) 1.35+0.03 


(Received June 10, 1943) 


The energies of some y-rays have been determined by the method of semicircular focusing 
of Compton recoils in a magnetic spectrograph. The results thus obtained are as follows: 


1.17+0.02, 2.65+0.06 0.94+0.02 


Ga™ we La“? 
2.04+0.04 


The two quanta emitted in the disintegration of Ga™ are present with equal intensity, 
which suggests that they may be in cascade. The y-ray activity of Ga™ was followed for 100 


INTRODUCTION 


HE energies of the y-rays emitted in the 
disintegration of several radio-elements 
have been measured by means of a magnetic 
spectrograph which has been previously de- 
scribed.! Compton recoil electrons are focused by 
a magnetic field, and coincidences are then 
observed as a function of Hp. The radius of 
curvature of the path of the recoils is 5.50 cm. 
Because of the fact that slow electrons are 
heavily absorbed and scattered from the focused 
beam by the walls of the counters and by the 
argon-alcohol counter mixture which is present 
throughout the magnet box, it has been found 
~ *At present at the Radiation Laborat Massachu- 
setts Institute of Technology, Cambridge, assachusetts. 


1C. E. Mandeville, Phys. Rev. 62, 309 (1942); Phys. 
Rev. 63, 387 (1943). 


hours and was found to decay with a half-period of 14.25+-0.20 hr. 


advisable to employ double coincidence counting, 
using the counters 7; and 7: (Fig. 1, reference 1), 
when making observations on quanta of an 
energy less than 1 Mev. Triple coincidence 
counting, with the counters 7), T2, and 73, is 
especially suitable for obtaining end points of 
distributions in regions of higher energy (greater 
than 1 Mev), since the gamma-ray background 
is then very small. The absorption of the slow 
electrons also leads to a lower limit of satis- 
factory measurement. The recoils arising from 
y-rays of energy less than about 0.5 Mev are 
not observable. 

The treatment of the background count and 
the corrections which are applied to data ob- 
tained with this spectrograph have been previ- 
ously outlined' and are obviously the same, 
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° 2 4 6 
Hp x10” OERSTED- CM 
Fic. 1. Momentum distribution of the Compton recoils 
of the y-rays from Sc**. The Sc* of this curve was produced 
by a bombardment of CaF; by deuterons of energy 12 Mev. 
riple coincidences in the counters 7;, 72, and 73 were 


recorded 


whether double or triple coincidences are 
recorded. 

In all instances under discussion in this paper, 
observations were begun at a time not less than 
twenty-four hours after cessation of bombard- 
ment. 


Sc** 


The y-rays from Sc** have been discussed in a 
preliminary report.? The momentum distribution 
of the Compton recoils of the y-rays from Sc*® 
produced by a bombardment of calcium fluoride 
by deuterons of energy 12 Mev is given in 
Fig. 1. The energy of the y-ray taken from the 
end point of the distribution is found to be 
1.35+0.03 Mev. The ordinates of the distribu- 
tion were found to decay with the established 
44-hr. half-period* of Sc**. In obtaining the curve 
of Fig. 1, triple coincidences were recorded. The 
electrons of low energy are seen to be heavily 
absorbed on the interval 0 <Hp <2000. 

Pure titanium dioxide was bombarded by fast 
neutrons from the reaction .Be-d—n. The mo- 
mentum distribution of the Compton recoils of 
the y-rays from Sc** formed by the reaction 
Ti**(n,p)Sc** is given in Fig. 2. The quantum 
energy obtained from the end point of the 


*C. E. Mandeville, Phys. Rev. 62, 555 (1942). 
*H. Walke, Phys. Rev. 57, 163 (1940). 


C. E. MANDEVILLE 


Fic. 2. Momentum distribution of the Compton recoils 
of the y-rays from Sc**. The Sc* of this curve was produced 
by the reaction Ti**(,p)Sc**. Double coincidences in the 
counters JT; and 7: were ed. 


distribution is again found to be 1.35+0.03 Mey, 
Double coincidences were recorded. 

Hibdon, Pool, and Kurbatov‘ have likewise 
found Sc** to be a gamma-ray emitter, and their 
absorption experiments indicate a quantum 
energy of 1.33 Mev. 


Ga” 


Pure metallic gallium has been irradiated by 
slow neutrons. The y-ray activity was followed 


ed 


ACTIVITY — ov,/ SEC. 


° 20 4 60 60 100 
"HOURS 


Fic. 3. Decay of gamma-ray activity of Ga™, The 
“eer taken from the slope of the curve is 14.25+0.20 
r. 


*C. T. Hibdon, M. L. Pool, and J. D. Kurbatov, Phys. 
Rev. 63, 462 (1943). 
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for one hundred hours with a Lauritsen electro- 
scope, and the half-period taken from the slope 
of the decay curve of Fig. 3 was found to be 
14.25+0.20 hours, in good agreement with the 
value 14.1+0.2 hr. obtained by Sagane.* The 
momentum distribution of the Compton recoils 
of the y-rays from Ga” is given in Fig. 4, where 
it is seen that two quanta of energies 1.17+0.02 
Mev and 2.65+0.06 Mev are emitted with equal 
intensity, suggesting that they may be in cascade. 
It is of interest to note that Sagane et al.,® 
using cloud chamber and magnetic field, have 
found the K-U extrapolated limit of the 6-rays 
from the negatron emitting Ga” to be 1.71 Mev, 
whereas absorption experiments by Livingood 
and Seaborg’ have indicated a maximum #-ray 
energy of 2.6 Mev. The disagreement between 
the two results might be related to secondary 
effects arising from the 2.65 Mev y-ray from 
Ga”. The influence of Compton recoils from 
y-rays of high energy upon absorption measure- 
ments of the energies of primary 8-rays has been 
discussed in several recent publications.*" 


2 4 6 8 


Hp x 10” OERSTED-cM 


Fic. 4. Momentum distribution of the Compton recoils 
of the y-rays from Ga”. The two electron groups corre- 
spond to two equally intense quanta having energies of 
1.17+0.02 Mev and 2.65 40.06 Mev. 


*R. Sagane, go Rev. 55, 31 (1939). 
*R. Sagane, S. Kojima, and G. Miyamoto, Proc. Phys. 
dG. T. Scaborg Rev. Mod. Ph 
. J. Livin and G. T. 
12, 30 (1940) 


N 
He 


° 3 4 
Hp x 10° OERSTED- cM 


Fic. 5. Momentum distribution of the Compton recoils 
of the y-rays from W"*", 


The radioactive tungsten isotope of half-period 
24.1 hours" was produced when pure powdered 
tungsten was irradiated by slow neutrons. The 
momentum distribution of the Compton recoils 
of the y-rays from W'* is given in Fig. 5. The 
quantum energy as determined from the end 
point of the distribution was found to be 0.94 
+0.02 Mev. The energy of this-y-ray has been 
previously reported to be 0.87+0.03 Mev by 
lead absorption," and aluminum absorption of 
the Compton recoils of this y-ray has led to a 
quantum energy of 0.90 Mev.” From the shape 
of the curve of Fig. 5, it is clear that the y-ray 
of energy 0.94 Mev is the only y-ray present 
with an energy greater than 0.5 Mev. 

Evidence for the emission by W'*’ of y-rays 
of an energy less than 0.5 Mev has also been 
obtained." 


The y-rays emitted by the 40-hr.“ La" have 
been discussed in a previous report.'! The y-ray 
activity of pure La,O; irradiated by slow neu- 
trons was followed for 160 hours and was found 
to decay with a half-period of 40.8+1.0 hr. The 
decay curve is given in Fig. 6. The momentum 


" Kasimir Fajans and W. H. Sullivan, Phys, Rev. 58, 
276 (1940). 

#2 Arnold F. Clark, Phys. Rev. 61, 242 (1942). 

3G. E. Valley, Phys. oy 59, 686 (1941). 

“4K. E. Weimer, M. L. Pool, and J. D. Kurbatov, 
Phys. Rev. 63, 67 (1943). 
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*M. Deutsch, Phys. Rev. 61, 672 (1941). 
*M. Deutsch, J. R. Downing, L. G. Elliot, J. W. 
Phys. Irvine, Jr., and A. Roberts, Phys. Rev. 62, 3 (1942). 
Lise Meitner, Phys. Rev. 63, 73 (1943). 
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ACTiviTy- oiv/sec. 


0 2 40 #60 60 100 120 140 160 160 


HOURS 
Fic. 6. Decay of -ray activity of La,O; irradiated 


> | slow neutrons. The half-period taken from the slope 
the curve is 40.8+1.0 hr. 


distribution of the Compton recoils of the y-rays 
from La"® is to be seen in Fig. 7, and the quantum 
energy of the single y-ray present was found to 
be 2.04+0.04 Mev, in good agreement with the 
previously determined absorption value of 2.00 
+0.05 Mev." 

It is not as yet certain whether La may 
be regarded as a source of absolutely mono- 
chromatic 2.04-Mev radiation, since recent coin- 
cidence experiments!® by W. C. Peacock of The 
Massachusetts Institute of Technology seem to 
indicate that a y-ray of energy less than 0.3 Mev 
is in cascade with the 2.04-Mev quantum. These 
coincidence experiments were only of a prelimi- 
nary nature, however, and are not yet to be 
regarded as conclusive. 


OTHER REMARKS 


Sc*® and La"® appear to be, for the most ‘part, 
emitters of monochromatic gamma-radiation. It 
is suggested, therefore, that the y-rays of these 
elements might be used for the determination 
of the exact shape of the momentum distribution 
of the Compton recoils of a single y-ray at the 
outset of experiments in which spectrographs of 


18M. Deutsch, private communication. 


He 10° OERSTED-cm 


Fic. 7. Momentum distribution of the Compton recoils 
of the y-rays from La™®, 


this type or related types are employed. A 
knowledge of the shape of the distribution 
arising from a single y-ray is particularly 
desirable, since a misunderstanding of the nature 


of the distribution may lead to an erroneous 
interpretation of results obtained from an exami- 
nation of emitters of complex gamma-ray 


spectra.' 
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The Zeeman effect of the rare gases is discussed from the standpoint of the coupling scheme 
of Shortley and Fried. The matrices for this type of interaction are calculated and also the g 
factors from the transformations to LS coupling for the p*f configurations. These enable us 
to determine the Paschen-Back patterns to be expected. The calculated patterns are compared 
with observed structures and the agreement is very satisfactory in most cases. 


E have become so accustomed to hearing 

the phrase “transition from LS to jj 
coupling” that whenever we think of the varia- 
tion in coupling in a series of levels, we auto- 
matically think of the change taking place from 
the LS condition toward a jj coupling limit. 
' The spectra of the rare gases afford examples 
which show that this simple picture is inade- 
quate. The work of Shortley and Fried' has 
shown that in the spectra of Cu II and A, a very 
striking example of a special type of intermediate 
coupling exists. This may be shown by calcula- 
tions which neglect all parameters in the energy 
matrix excepting F, and ¢,. This is equivalent 
to assigning a J value to the core (p*) and using 
only one electrostatic parameter in addition. 
The result is to break up the spectrum into two 
widely separated groups (?P3;2 and *Py2 of the 
p® core) each of these in turn being broken up 
into a group of closely-packed levels. These 
closely-packed levels in turn exhibit the pe- 
culiarity of being levels with multiple J values. 
Shortley and Fried! applied their results only to 
the spectrum of argon among the rare gases, 
but since then a considerable amount of data on 
the Zeeman effect has been published which 

helps to substantiate this picture. 

Some years ago,” the author had occasion to 
investigate the genealogy of the p°p configura- 
tions of the rare gases as it could be determined 
from the g factors. Satisfactory curves could be 
drawn for the values known up to that time 
showing the transition from LS coupling toward 
jj coupling, but the slopes of some of the curves 
seemed to be too abrupt at the jj end. Since 


1G. H. Shortley and B. Fried, Phys. Rev. 54, 749 (1938). 
? J. B. Green, Phys. Rev. 52, 736 (1937). 


that time, a large amount of data on the p*d 
and some data on the p'f configurations have 
been given.’ It was primarily the data on the 

pf configuration which led to the considerations 
discussed in this paper. 

Table I gives the matrices for the transforma- 
tion from LS to jj coupling of the pf configura- 
tion, together with the g values to be derived 
from it by performing the operation 


g(aJ) = a(SLJ) (ySLJ|aJ)?. (1) 


Thus, to find the g value corresponding to a 
particular row JJ’, square the numbers appear- 
ing in that row and multiply each of the numbers 
thus obtained by the g(SLJ) occurring below 
the LS designation at the head of the column in 
which the number lies. Add all the terms thus 
obtained and this will be the g value for jj 
coupling. This same method applied to the 
transformation matrix between any other type 
of coupling and LS coupling will yield the g 
value for that type of coupling. 

It is, of course, not necessary to perform these 
laborious calculations for the sole ‘purpose of 
determining the g factors in any limiting case of 
coupling. Closed formulas for the calculation of 
g factors have been given for a number of these.‘ 
But for the calculations of intensities and energy 
levels in intermediate fields, the complete matrix 
must be known.$ 


(a) Neon: . Green nag: G . Peoples, Phys. Rev. 54, 
602 (1988). ® Argon: Green and B. ried, Phys. 
Rev. 54, 876 (1938). op rg Green, Bowman, and 
Hurlburt, Phys. Rev. 58, 1 (1940). (d) Xenon: Green, 
Hurlburt, and Bowman, Phys. Rev. 59, 72 (1941). 

* See, e.g., Pauling and Goudsmit, The Structure of Line 

icGraw- ill, New York, 1930). 
G. , Phys. "Rev. 61, 537 (1942) has used the 
engine method to calculate g formulas for this 
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TABLE II. Matrices of electrostatic energy F; in jj 


TABLE I. Matrices for transformation from LS to jj 
coupling for 


coupling for 
— Fo+ 


1.200 


4c 


2(15)' _10(35)! 


v3 —6 


—(21)§ 2/7 


NIK NIW 


NIN 


1.250 1.000 


7 
6/7 —12./5 6(21)8 

9(15)! 


F; 


5(35) —2(105)! 


NIK NIK NIW 
NIG NIN Nin 


—3/3 5(21)! 8(15)§ 677 


36 —4/7 275 4(21)! 


a 


750 1.083 1.333 1.000 . F; 


24(14)! Fs 


Nin 


2/7 —(14)' (2199 ~12F; 


a 
Nie 


0.667 1.167 1.000 
TABLE III. Transformation matrices from jj to 


2D, gif) intermediate coupling. 


1 0.500 
3a’ 


1 
(21) 


g(LS) 0.500 


Table II gives the matrices of the electrostatic 
parameter F:, the only one to be considered in 


ile of coupling, but they give the same result as the 


151)l2]s2} formulas that are given in Pauling and 
dsmit (reference 4), or Candler, Atomic Spectra 
(Cambridge Univ. Press, 1937, p. 143-145). 
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this example, and if the determinants of these 
matrices are solved, they yield the following 
results, after adding the electromagnetic param- 
eter ¢, to the higher set (c and d), and subtracting 
it, from the lower set (a and 6). 


Upper set 4c, 3c, 3d, 2d= — Fyt+f>. 


Lower set 5a, 4b’= —Fo—}f,— (1) 
4a’, 3a’ = — Fyo—3f,+10F:, 
3b’, 2a’ = — 3 Fo, 


2b’, 1b = — Fo—}tp—12F >. 


These solutions, inserted in Table II, give the 
transformation matrices for jj coupling to this 
particular variety of intermediate coupling and 
are given in Table III. When we apply these 
transformations to the matrices of Table I we 
arrive at Table IV, which gives the transforma- 
tion from ZS to this intermediate type of 
coupling. Listed with Table IV are the g values 
for this special coupling scheme. A striking 
difference is shown between the g values for the 
intermediate case and those for jj coupling, even 
though they seem to be quite close together at 
first glance. Perhaps even more striking is the 
difference of phase between several homologous 
elements of the two sets of matrices. 

A study of the p'f configurations of the rare 
gases from the standpoint of Eqs. (1) shows some 
interesting points. 


Strangely enough the spectrum of neon shows 
the least conformity with the above results. The 
levels labeled W (except 4W) are separated from 
the others of the lower (?P3;2) set by about 100 
cm, Except for 4W, all the members of the 
upper set are piled together, as indicated by (1). 
The author is of the opinion that interaction 
with the 5p configurations is not sufficiently 
large to account for such large perturbations, 
and suggests that some revision of these assign- 
ments is necessary. 


Argon 
As indicated by Shortley and Fried’ the pf 
configurations of argon represent an almost 


*W. F. M . Humph Bur. Stand. J 
10, 427 ( 540). C and 
Bur. Stand. J. Research 10, 139 (1933) 
521). J. Humphreys, J. Nat. Bur. Stand. 
20, 17 (1938) (RP 1061). 
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perfect example of the application of Eqs. (1). 
Each of the levels of the lower (?P 3,2) set can be 
assigned a double J value, and the spacings of 
the four observed levels are in almost perfect 
accord with the theory. The upper (?P,2) set 
should consist of four levels lumped into one, 
but evidence of splitting into two levels is 
noticed. 


Krypton 


In krypton, we see for the first time the 
splitting of the levels of the lower (*P3,2) set into 
pairs. This appears as a separation of the lowest 
level (J=1, 2) into two separate levels (X, Z); 
and of the third (J=2, 3) into two separate (but 


TABLE IV. Matrices for transformation from LS to 
intermediate coupling. 


gj) 
sa’ | 1 | 1.200 J=s 
g(LS) 1.200 

%G, 1G, g(Int.) 
| VS 373 2 1.195 
4c’ | -3 2/3 1,083 


g(ZS) 1.050 1.250 1.000 


OF, g(Int.) 
3a’ | 06/21 0.964 
36’ | 0 1.175 
xt | sy7 ary | 42 1.206 
3a’ | 213 -3(21)) 0 0.821 
g(LS) .750 1.083 1.333 1.000 

Ds g(Int.) 
| -272 273 0.844 
2a’ | 2/5 —(10)) (15)! 0.889 
g(LS) .667 1.167 1.000 

g(Int.) 
1’ 1 | 0.s00 J=1 


J. B. GREEN 


Fic. 1. (a) \\ 6178-9, Xe 1s2—SY, 5X. (b) AX. 7584, Xe 3d,’—5ST, 5Z. 


(The lines to the right of the figure are Kr.) (c) AA 7472-4, Xe 3ds—SY, 5X. 
Above are shown microphotograms of original plates; below are shown 


calculated patterns. Full lines are 
parallel 


very close) levels (Y, 7). 5X and 5Z are per- 
turbed by the upper terms of p°3p and are 
therefore more widely separated than we should 
expect. No levels of the upper (?P1/2) set are 
known for krypton. 


Xenon 


This is probably the best known spectrum as 
far as the p'f configurations are concerned. The 
members of the lower set are all separated with 
the exception of the highest (W), but lines 
involving these levels show Paschen-Back effect, 
indicating that the W levels are complex. All of 
the levels of the upper (?P1/2) set lie in the 
negative energy region and none has been found. 

Lines involving the p'f configurations in neon 
and argon did not appear on the plates. A few 
of these lines appeared on the krypton plates, 
and many of them on the xenon plates. 

Examples of the Paschen-Back effect of xenon 
chosen for this paper were the p*5f(2b’, 1b) and 
the p*5f(5a, 4b’) known as SY, 5X and 57, 5Z, 
respectively. These particular levels were chosen 
because each pair contains one pure level 


dicular polarizations ; dotted lines 


tions; dot and dash, mixed polarizations. 


(SX =*D, and 5T7=*F;). The results are shown 
in the accompanying figure. The method of 
making these calculations has been outlined ina 


TABLE V. g values for p*p of rare gases. 


3 

& 


J=1 pro Ne 


1.984 
A 1.985 1.90 JJ 1.500 
Kr 1.898 1.834 1.795 1.795 = Int. 1.778 
Xe 1.852 1.728* 1.801 
J=1 pi Ne 0.669 0.974 
A 0.838 1.01 JJ 1,333 
Kr 1.004 1.034 1.041 1.014 Int. 1.056 
Xe 1.022 0.903* 1.036 
J=1 Ps Ne 0.999 0.685 
pe A 0.819 0.61 JJ 0.667 
pe Kr 0.647 0.648 Int. 0.611 
Xe 0.790* . 
J=1 pe Ne 1.340 1.397 
pe A 1.380 1.45 JJ 1,500 
Da 1.452 1.401 Int. 1.556 
Xe  =1.552* 
J=2 pa Ne 1.301 1.184 
Ds A 1.260 JJ 1067 
pe Kr 1.181 1.158 Int. 1,067 
ps Xe 1.195 
J=2 pe Ne 1.229 1.360 
pe A 1.305 1.42 JJ 1333 
Pe Kr 1.388 1.403 1.403 1.411 Int. 1.433 
fe Xe 1.379 1.347 1.395 1.386 
J=2 Ds Ne 1.137 1.112 
ps A 1.112 1.09 JJ 1.167 
pr Kr 1.099 1.107 Int. 1.067 
pe Xe ——:11.106 1.123 1.103 


* These levels from different configurations are close enough to 
perturb each other. 


J: 


Je 
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TABLE VII. g values for p*f configurations of the rare gases. 


Ne 0.752 0.797 JJ 0.833 

J=t A 0.877 0.768 int. 0.833 
jai ds Ne 0.860 0.812 

A. 0.768 0.813 (1.186)? JI 1110 

Int. 0.833 


A 400 JJ 1.333 
Kr 1.098 1.348 1.355 1 Int. 1.067 
Xe 

J=2 Ne 
A 0.908 0.941 1.107 JJ 1,121 
Kr 1.006 0.965 0.954 1.005 Int. 0.978 
Xe 1.073 0.987 0.980 

J=2 ds Ne 1.356 1.322 1.298 
A 1.437 1.387 1.206 JJ 1.067 
Kr 1.295 1.318 1.315 Int. 1 
Xe 1.376 1.196? 1.303 1.298 

J=2 si” Ne 1.242 1.230 1.251 
A 0.987 1.265 1.264 JJ 1.289 
Kr 0.899 Int. 1.300 

J=2 si” Ne 0.781 0.783 
A 1.057 0.802 0.777 JJ 0.767 
Kr 1.169 Int. 0.756 
Xe 1.274* 

J=3 Ne 1.034 1.040 1.093 

- A 1.077 1.076 1.052 JJ 1.068 

Kr 1.050 1.073 1.094 1.037 Int. 1.036 
Xe 1.026 1.076 1.078 1.081 1.082 

J=3 dy’ Ne 1.249 1.248 
A 1.199 1.245 JJ 1.237 
Kr 1.243 1.254 1.231 1.227 Int. 1.270 
Xe 1.263 1.246 1.225 

J=3 si” Ne 1.125 1.116 
A 1.133 1.127 1.098 JJ 
Kr 1.140 Int. 1.111 
Xe 1.126 

* Should probably be labelled si’. 


previous paper.’ Excellent agreement is seen to 
exist between calculation and experiment. In 
the case of 5Y, 5X the.theoretical g values are 
1.10 and 0.50. The values calculated from the 
pattern are 1.11 and 0.50. The difference between 


7]. B. Green, Phys. Rev. 59, 69 (1941). 


JJ Int. 
J=1 x Kr 0.50 0.61* 0.500 0.500 
J=1 xX Xe 0.504 0.50 0.50 0.500 0.500 
J =2 V Xe 0.86 0.87 0.87 1.047 0.844 
J =2 Y Xe 1.11 1.10 1.09 0.897 1.100 
J =3 U Xe 1.18 1.17 0.997 1.175 
J=4 Ws Xe 1.19 1.181 1.195 
J=4 Zz Xe 1.02 1.034 1.022 


* Perturbed by 4p*6p. 


the observed and theoretical g value of the 5Y 
level can be accounted for by the fact that the 
5Y and 5X levels are separated by 3.6 cm in- 
stead of falling on top of each other. (See Fig. 1.) 

With respect to the V, U (30’, 2a’) interaction, 
a serious discrepancy must be noted. Although 
the g values fit the theoretical calculations 
extremely well, the magnetic levels should not 
perturb each other very much because the 
interaction to be expected in this case almost 
cancels out. Yet the observed patterns are ma- 
terially perturbed, and enough so to be somewhat 
greater than experimental error. Either the 
theory described here is too simple to be used in 
this particular case, or some perturbations from 
outside configurations are at work. 

Table V gives a summary of the available 
data on the g factors of the rare gases for p*p. 

The tendency toward the limit for intermediate 
coupling is particularly evident for pio. 

Table VI gives a summary of the data for p°d 
configurations. 

Table VII gives the summary of experimental 
data for p*f configurations. 

It is in this configuration that the evidence in 
favor of the particular type of coupling assumed 
in this paper is most striking. 
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Owing to cumulative ionization, in calculating the bal- 
ance between the impact ionization and the recombination 
in a homogeneous electric field, three possible equilibrium 
values for the number of free electrons are obtained. It 
results in an S-like current vs. field strength curve. The 
breakdown is interpreted as an abrupt transition from the 
lower branch of the curve to the upper one. An expression 


for the breakdown field strength is obtained. It turns out 
to be of the correct order of magnitude. In a certain region 
of current densities, no stable homogeneous solutions exist, 
and the current is concentrated in the shape of a narrow 
discharge canal—a positive column. A non-linear differen. 
tial equation for the distribution of free electrons in its 
cross section is obtained. 


1. 


N considering the electric current through a 
gas in a homogeneous field, at not very low 
pressure and at sufficiently great distance be- 
tween the electrodes, the boundary regions and 
those in the vicinity of the electrodes may be 
neglected. The number of free electrons obtained 
as a result of impact ionization is determined in 
this case by the equilibrium between ionization 
and recombination at every point: 


(1) 


Here n stands for the estimated number of free 
electrons per unit volume; JN is the number of gas 
molecules; a@ is the ionization coefficient; p—that 
of recombination. The number of positive ions is 
supposed not to differ essentially from the num- 
ber of free electrons. 

- At any given field a and p have definite values, 
and Eq. (1) gives a definite value of which 
increases with increasing field strength. Because 
of this increase in m, the current also increases 
continuously. 

In reality, however, a continuous increase of 
the current proceeds only up to a certain value of 
the field, while at the further increase of the field, 
an abrupt change of the current takes place, 
that is, breakdown or the ignition of the gas dis- 
charge. Such an abrupt change may be obtained, 
if the cumulative ionization, that is the impact 
ionization of molecules, previously excited by the 
electron impact is taken into account in Eq. (1). 

For the sake of simplicity we shall draw our 
attention only to two excited levels (1 and 2) 
between the ground level of the electron in a 
molecule (0) and the continuous spectrum. Let us 


an N =pn?. 


denote the number of excited molecules per unit 
volume by m, and mz, respectively. Usually 
nXKn2“Kn, and in the approximate evaluation of 
n, the transitions to higher levels may be neg. 
lected. This gives 


Bown N =(y1 (2a) 


Here 80; is the probability of collision of second 
kind between the excited molecule and the free 


electron. 
In a similar way for the second level, if we 


neglect the ionization, 
(Bo2N +B12m1)n = (y2 m2, 


where f12 is the probability of secondary excita- 
tion 1—2, all the other designations are similar to 
those of Eq. (2a). Finally, for the number of 
collisions leading to ionization: 


(2b) 


(aoN =pn’. (2c) 


AE 
Fic. 1. Characteristic curve. 
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Here a; is the probability of the ionization of the 
molecule in the ith excited state. Equations (2) 


give for m: 
(3) 
where 
a=pb1b2, Ci=pryiv2; (4) 
= (205152 + 8010152 1202) N, 
JN, (5) 
d=aoyi172- 


The number of free electrons capable of 
ionizing and exciting increases rapidly with the 
increase of strength of the applied field E. In a 
simplest case, when the mean free path of the 
electrons does not depend on their velocity, this 
number is proportional to! 


ox) 


where ¢; is the energy of excitation or of ioniza- 
tion, respectively, / is the mean free path of the 
electrons; C1. The coefficients a; and 8% must 
depend on the field strength approximately in 
the same way, and in comparison to the sharp 
exponential dependence, y; and 6; may be con- 
sidered independent of the field. Consequently 
a, b;, and c,; also do not depend on the field, 
while bz, c2, and d in weak fields are very small 
(be&b; and cc) and increase strongly with the 
increase of the field strength. 

The shape of the curve n=n(E) according to 
Eq. (3) depends on the behavior of the coeffi- 


f~exp ( 


cients of this equation. If ce, increasing with the * 


increase of the field strength, reaches the value of 
¢; before reaches the value of then a 
monotonously increasing curve results. If, how- 
ever, b2.=b, already at ¢2<c,, then an S-like 
curve is obtained (Fig. 1). 

In reality, because of the fact that the distri- 
bution of the velocities of electrons in a homo- 
geneous field decreases at a greater rate than the 
Maxwell distribution, the second case takes 
place. Thus, with the dependence prescribed by 


1M. J. Druyvesteyn, Physica 10, 61 (1930); 1 — 
(1934); 3, 65 (1936). B. Davydov, Physik. Zeits. Sow} 
union 8, 59 (1935); 9, 433 (1936). P. M. Morse, W. P. Allis, 
and E. S. Lamar, Phys. Rev. 48, 412 (1935). 


ELECTRIC BREAKDOWN 


Eq. (6) we get 
(5’) 


In cz we have neglected Bo2a271 in comparison 
with Boary2. Usually €9:>¢€, and, due to this the 
exponential factor in the first of these quantities, 
is smaller than in the second. Now, by the order 


of magnitude 
but 


€12€2 
(2c —); 
consequently at 
_ The solution of the Eq. (3) may be approxi- 
mately substituted for by the solution of the 
square equations. At small m, neglecting the term 
with n', we get 


1 
"a 
(7) 


The lower sign before the square root gives the 
first lower branch of the curve, the upper, the 
second branch. Both of the branches join at a 
point A (Fig. 1), when the expression under the 
square root disappears, and this gives the lower 
turn of the S-like characteristic. 

‘At large values of m, d in Eq. (3) may be 
neglected. Then 


1 


For the middle of the second branch, Eqs. (7) 
and (8) give approximately the same value: 


n= 


If we increase the applied field strength, 
moving along the lower branch of the curve, then 
at a point A or close to it, an abrupt transition to 
the upper third branch—a breakdown or an 
ignition—should take place. We get the condi- 
tion for the breakdown potential: 


(9) 


or approximately 
C=C). (9’) 


According to Eqs. (4) and (5’) this gives 


PY1- (10) 
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With known approximate values of the quanti- — 


ties involved, this gives the correct order of 
magnitude of the breakdown field. 

If we take into account not only the two ex- 
cited levels, as was done above, but a greater 
number of them, then in the same way, instead of 
Eq. (3), an equation of a higher degree should be 
obtained. However, for small values of m and for 
the region close to the lower turn of the charac- 
teristic, we may also in this equation neglect the 
higher powers of m, as we neglect the term with n® 
in solving Eq. (3). Because of this, the curve 
n(E) at small values of m does not change es- 
sentially. Conditions (9) and (10) for the break- 
down field shall particularly remain valid. By the 
way, even the second level does not enter 
Eq. (10), despite the fact that it is necessary in 
obtaining the turn of the curve. 

Equation (10) is in an approximate agreement 
with the Paschen’s law. Indeed, all the coeffi- 
cients a; and By depend only on the product El 
[compare Eq. (6) ], i.e., only on pd, where p is the 
gas pressure, d is the distance between the 
electrodes. The fact that » appears in Eq. (10) 
also through N cannot have any noticeable effect, 
since the dependence of p and a; on Bix is much 
stronger. 


The second middle branch of the characteristic 
is not stable. In Fig. 1, to the right of the curve, 
i.e., above the second branch, ionization over- 
weighs the recombination. To the left, i.e., below 
the second branch, recombination overweighs. 
Because of this, close to the second branch, 
random deviations of m from the value prescribed 
by Eq. (3) shall not only not disappear, but, on 
the opposite, shall increase with time, Having 
included in series with a gas gap a very large 
constant resistance, it is always possible to fix not 
the field strength but practically the current 
intensity. The instability of the second branch 
of the characteristic leads to the fact that in the 
region of currents corresponding to m lying 
within the range from point A to point B no 
stable homogeneous states exist. Random fluctu- 
ations lead to the concentration of the current in 
the form of a narrow discharge canal—a positive 
column—in full accordance with the phenomena 


observed. 


(an?+bn -— 


B. DAVYDOV 


The existing theory of the positive column? 
does not take into consideration the cumulative 
ionization. It may be taken into account, if 
Eq. (3) is complemented by the diffusion terms, 
Because of the high concentration of the free 
electrons in the positive column, there exists a 
Maxwellian distribution of their velocities with 
an abnormally high temperature of the electrons, 
which depends on the longitudinal field. In this 
case, in the radial direction a quasi-neutral diffy- 
sion of electrons and positive ions with the 
ambipolar coefficients of the diffusion 


(11) 
e 

is obtained. Here uz and 7+ are the mobilities 

and the temperatures of the free electrons and 

the positive ions, respectively. 

Now in Eq. (3), to the decrease of the number 
of free electrons due to the recombination, we 
must add their decrease as a consequence of the 
ambipolar diffusion, i.e., we must make a 
substitution: 


This leads to the non-linear differential equation 


1 
det z dz 
=byn?+con+d. (12) 


Outside of the discharge canal, estimated solu- 
tions of this equation must coincide with m,, i.e, 
with the first root of Eq. (3) (Eq. (7) with the 
lower sign before the root sign): 

With 


2, 


On the other hand, the boundary condition at 
z=0 will be dn/dz=0. Here n(0) itself must lie 
between the second and the third root of Eq. (3). 

This sort of stable solutions of Eq. (12) must 
fill also the region of the current densities corre- 
sponding to the interval between the points 4 
and B. 

A more detailed account of this work will be 
published in the Journal of Physics U.S.S.R. 


*W. Schottky, Physik. Zeits. 25, 342 and 635 (1924); 
L. Tonks and I. Langmuir, Phys. Rev. 34, 876 (1929). 
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The Propagation of Order in Crystal Lattices* 


SEPTEMBER 1 AND 15, 1943 


J. AsHKIN AND W. E. Lams, Jr. 
Pupin Physics Laboratories, Columbia University, New York, New York 


Zernike’s problem of the propagation of order in a binary 
crystal alloy is discussed by means of the matrix formalism 
recently developed for treating cooperative phenomena. 
It is proved generally that the existence of long range order 
over any temperature range implies a degeneracy of the 
maximum characteristic value of the fundamental matrix 
over the same temperature range. For a two-dimensional 
crystal, the special form of the matrix theory developed by 
Kramers and Wannier has been used to obtain explicitly 
the probabilities for finding an A or B atom at any lattice 
site if it is known that there is an A (or B) atom at a certain 
site (correlation probabilities, or intermediate range order). 
The results are obtained in the form of power series valid 
at low temperatures by a perturbation treatment of the 
highest characteristic values and characteristic vectors. 
It is seen that the maximum characteristic value is doubly 
degenerate from T=0 to some finite temperature, consist- 
ent with the existence of long range order at low tempera- 


(Received June 10, 1943) 


ture. The absence of long range order at high temperatures 
is proved by showing that the maximum characteristic 
value is non-degenerate at sufficiently high temperatures. 
By comparing the solutions of Zernike’s approximate equa- 
tions for the correlation probabilities in two dimensions 
with our exact solution, and his expressions for the long 
range order and energy (short range order) in three di- 
mensions with those given by van der Waerden, we find 
that Zernike’s approximation is better for a two- than a 
three-dimensional crystal. The problem of propagation of 
order is generalized by an investigation of the correlation 
probabilities for more complicated configurations of fixed 
atoms. For example, the ordering influence of an adjacent 
pair of disordered atoms (“dipole”’) is found as a function 
of position in the lattice. It is found here, as might be 
expected, that the ordering influence falls off rapidly with 
the distance of the site from the dipole. Other cases are 
treated as well. 


N order to interpret the diffuse scattering of 
x-rays by a partially disordered crystal alloy 

it is necessary to know how the presence of one 
atom at a certain lattice point affects the proba- 
bility of finding the various kinds of atoms at the 
other lattice points. The older theories of the 
order-disorder problem, principally those of 
Bragg and Williams,’ and Bethe,’ gave only the 
long and short range order and did not consider 
the more general problem of finding the correla- 
tion proabilities for all distances, intermediate 
as well as long and short. Experimentally, several 
of these correlation probabilities were measured 
by Wilchinsky* for the alloy CusAu, while the 
problem of calculating them theoretically was 
first attacked by Zernike* in a paper entitled 
“The Propagation of Order in Cooperative 
Phenomena.”’ 

In Zernike’s work, the propagation of order is 
described by a partial finite difference equation 

* Submitted in partial fulfillment of the requirements for 
the degree of Doctor of Philosophy in the Faculty of Pure 
Science, Columbia University, New York, New York. 
Publication assisted by the Ernest Kempton Adams Fund 
for Physical Research of Columbia University. 

1W. L. Bragg and E. J. Williams Proc. Roy Soc. London 
145A, 699 (1934). 

*H. A. Bethe, Proc. Roy Soc. 150A, 552 (1935). 


*Z. W. Wilchinsky, Phys. Rev. 63, 223 (1943). 
‘F, Zernike, Physica 7, 565 (1940). 
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relating the probability of finding a certain kind 
of atom at a given lattice point with the prob- 
abilities for various atoms at the nearest neigh- 
boring lattice points. However, his equation is 
necessarily approximate because in order to 
obtain a recurrence relation, he has to treat the 
probabilities for the nearest neighboring atoms 
as at least formally independent of one another. 
Nevertheless, this method is so far the only really 
practicable one for the three-dimensional propa- 
gation problem and gives reasonable results for 
the thermal properties and the correlation. prob- 
abilities at low and high temperatures. One of 
the results of the present investigation will be to 
obtain a comparison in two dimensions between 
Zernike’s solution and the exact solution valid at 
low temperatures. In three dimensions, exact 
series expansions for the long range order and 
the energy at low temperatures have been given 
by van der Waerden® and have been compared 
with the series of Zernike. The comparison indi- 
cates that Zernike’s approximation is better in 
two dimensions than it is in three. In one dimen- 
sion Zernike’s solution can be shown to be exact. 
It therefore appears that the approximation 


5 v. d. Waerden, Zeits. f, Physik 118, 473 (1941). 
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becomes successively poorer as the dimension is 
increased. 

Although there was no further work directly 
following up that of Zernike, a new technique has 
been introduced into the study of the order- 
disorder phenomenon which shows some promise 
of providing a solution to the propagation 
problem. Recent work by Montroll,* Kramers 
and Wannier,’? Lassettre and Howe,® and 
Onsager,® has shown how the statistical treat- 
ment of binary alloys can in principle be reduced 
to the solution of the characteristic value 
problem of a certain matrix, the characteristic 
values being related to the partition function for 
the crystal and the characteristic vectors to the 
probability for finding various configurations of 
the atoms in the lattice. It is the object of this 
paper to show how these matrix methods can be 
used to give at least a partial solution to the 
propagation problem. 

The application of the matrix theory to real 
three-dimensional crystals has been held up un- 
fortunately by considerable mathematical dif- 
ficulties. The matrix whose characteristic value 
problem is to be solved is very difficult to handle 
even by perturbation methods, the essential 
complication arising from the fact that all its 
elements are different from zero and do not 
appear to vary in any regular manner from one 
position in the matrix to the next. These difficul- 
ties have so far compelled us to restrict our inves- 
tigations to two-dimensional crystals, where the 
general matrix theory has been developed by 
Kramers and Wannier in a more tractable form. 
Using a perturbation method, valid at low 
temperatures, we have obtained power series for 
the maximum characteristic value and for the 
corresponding characteristic vectors. By inserting 
these series into the appropriate formalism of the 
matrix theory we then discuss the propagation 
problem by giving the correlation probabilities 
as functions of temperature and position in the 
lattice, obtaining series for the long range order, 
the short range order, and the order of inter- 
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mediate range. The treatment is restricted to low 
temperatures, however, the perturbation theory 
at high temperatures being very complicated and 
difficult to carry through. 

The paper is divided into two parts. Part | 
gives a brief outline of the matrix theory in three 
dimensions and contains a simple demonstration 
of the fact that states of long range order are not 
possible over a range of temperatures unless the 
maximum characteristic value of the associated 
matrix is degenerate over this same range. Part 
II is devoted to the discussion of the propagation 
problem in two-dimensional crystals. 


I. LONG RANGE ORDER AND THE 
MATRIX METHOD 


1. Results of the General Matrix Theory in 
Three Dimensions 

In this section we shall collect the results of 
the general matrix theory which are applicable 
to the problem under consideration. The notation 
is identical with that used by Montroll.® 

Consider a binary substitutional alloy AB 
whose crystal form is simple cubic. Let the edges 
of the crystal be L, M, N in units of the lattice 
distance. The crystal can then be regarded as 
made up of L layers of MXN atoms each, with 
each layer capable of 2¥* distinct configurations, 
corresponding to the fact that each of the MXN 
sites in the layer can be occupied by an A ora 
B atom. With the symbol a; used for a typical 
one of the configurations available to the atoms 
in the jth layer, v(a;) is to represent the total 
potential energy of interaction between all the 
nearest neighboring atoms of the jth layer in 
configuration a;, and v(a;, aj41) is taken as the 
representative of the total energy of interaction 
between atoms in the jth layer in configuration 
a; with their nearest neighboring atoms in the 
(j+1)th layer in configuration a;,;. With this 
notation and with the periodic boundary condi- 
tions imposed by requiring layer 1 to interact 
with layer L,'® the partition function has the form, 


Z=2 exp — {v(a1) +0(az) 
aL 
+(a1, a2) +--+ +0(az-1, az) 
+v(az, a:)}/kT. (1) 


10 This restriction is not necessary, serving merely to 
simplify the final formulae. 
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Making the substitution | 
V (ai, = +0(a;, 
we find that (1.1) becomes 


L 
Z= > > Il exp (- V(ai, ai41)/RT) 9 (1.2) 
fai} fay) 


with the understanding that 
Writing 
Ma, a’) =exp [— V(a, a’)/kT] 


we define a symmetrical matrix J whose 
2™Nx2™N rows and columns are labeled by 
a, a’, in accordance with the possible 
configurations a. The characteristic values \, and 
characteristic vectors t, (r=1, 2, ---, 2”%) of 
this matrix J are fundamental for the theory. 
For, as we shall see, if we can solve the char- 
acteristic value problem 


(1.4) 


we shall have the answer to all questions of 
thermodynamic and statistical interest for the 
crystal, including the order-disorder phenomenon. 

Acharacteristic vector has components, 
¢,(a), one for each configuration a. If we assume 
the A, and ¢,(a@) are known, the matrix element 
(3) can be written as 


(1.3) 


exp [— V(a, a’)/kT]= Artie(a) (1.5) 
provided that the ¥, are normalized to unity; 


Le = brs. (1.6) 


Substituting (1.5) into (1.2) and applying (1.6) 
one can easily prove that 


/1.7) 


r=1 

Thus the evaluation of Z, and consequently 
the investigation of thermodynamic quantities 
like entropy, energy, and specific heat, are 
reduced to the characteristic value problem (1.4). 
Since the number of layers L is very large for an 
actual crystal, the problem (1.4) can be restricted 
to that of obtaining the highest characteristic 
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value Amax and its corresponding characteristic 
vectors. If Amax is d-fold degenerate" (1.7) be- 
comes” 


(1.8) 


In the work of Montroll, and Kramers and 
Wannier it was assumed that Amax was non- 
degenerate for all temperatures with the ex- 
ception of possible isolated points. It was first 
pointed out by Lassettre and Howe?’ in their 
papers on binary solid solutions, that the maxi- 
mum characteristic value is actually doubly 
degenerate for a range of temperatures beginning 
at T=0, and they were able to demonstrate that 
within this range, a separation into two phases 
takes place. In the next section we shall show 
directly that this degeneracy of Amax is a neces- 
sary condition for the existence of long range 
order. 


2. Long Range Order and the Degeneracy of 2....x 


Suppose that a state of long range order exists 
in the AB alloy and that it is known that the 
atoms in layer 1 have a certain configuration a. 
Then the probability for finding the atoms of 
layer n, very far from layer 1, in a certain con- 
figuration a,, should depend explicitly on the 
configuration a. We shall show that if Amax is 
non-degenerate, the probability for the con- 
figuration a, is entirely independent of what 
configuration exists in layer 1, so that the non- 
degeneracy of Amax Will imply the absence of long 
range order. 

The probability P(a:,a,) of finding simul- 
taneously the configuration a in layer 1 and a, 
in layer m, is given by an expression similar to 
(1.2) where, however, the configurations a, and 
a, of layers 1 and m are held fixed. Using (1.5), 


" Although Amax itself may not be e, there may 
ual to Amax- As is pro or a special case in part II, 
will be even under such circumstances. A formula identical 
with (1.8) will then apply with d having an obvious mean- 
ing. In the future a statement that Amax is degenerate will 

be understood to include the possibility mentioned here. 
% It must be assumed here that the sum of the Lth 
powers of the canes characteristic values is negligible 
compared with Amax. Thus even if the ratios of the other 
characteristic values to Amax are small, say less than 1—x 
with x close to unity, there are ser 2% of them 
and there is no assurance that 2“* (1—x)*” is small com- 


pared to unity. 
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(1.6), and (1.7) we find 


oN 


s=1 t=1 


X teen) (2.1) 


To discuss the long range order, m and L—n 
must be taken as very large integers, the latter 
condition being necessary because of the bound- 
ary condition that layers 1 and L interact. This 
permits us to neglect all the characteristic 
values with the exception of Amax, and if the 
latter is assumed to be non-degenerate, (2.1) 


reduces to 
P(as, Gn) Umax (a1)? Umax(an)?. 


If P(ai, a.) is summed over all possible con- 
figurations a,, we obtain the probability P(a:) 
for finding layer 1 in configuration a. Since Umax 
is normalized to unity, this summation gives 
P (a1) =ttmax(a1)*. Similarly, P(an) = Umax (an). 
Thus P(a1, an) =P(a:)-P(a,). The probability 
for finding the configuration a, when it is known 
that the atoms of layer 1 are in configuration a, 
is equal to P(ai,a,) divided by P(a;). This 
quotient is completely independent of a:, being 
just equal to the probability for finding a, when 
no restrictions at all are put on the other layers. 
There is therefore no long range order under 
these circumstances. 

A further remark" concerning the interpreta- 
tion of this result must be made, in the light of 
a theorem of Frobenius to the effect that if all 
the elements of a matrix are positive, the maxi- 
mum characteristic value is likewise positive and 
in addition is non-degenerate. The elements of 
the matrix I, namely, exp [— V(a, a’)/kT], are 
all greater than zero and will evidently remain 
bounded away from zero so long as M and N 
remain finite. Consequently, for any crystal in 


P(a1, on) = 


(2.2) 


which the dimension L is very large compared 


with the dimensions M and N, vind will exceed 
by far the Lth power of any other characteristic 
value and states of long range order will not exist. 
On the other hand, if at least one of M and N is 
very large, and comparable to Z in magnitude, 
it may happen that some of the matrix elements 


% E. W. Montroll, J. Chem. Phys. 9, 711 (1941). 
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of I are close enough to zero to cause a prac. 
tical, if not actual, degeneracy in Amax; that is, 
there may exist other characteristic values }, 
whose Lth powers are not appreciably less than 


e. If Amax is d-fold degenerate, (2.1) becomes 


s=1 t=1 


P(ai, Qn) > Umax, 8(@1) Umax, 


x Umax, t(@n) Umax, ¢(a1) ’ (2.3) 


where the tmax,. are the characteristic vectors 
corresponding to Amax- Long range order will then 
in general be possible, since P(a:, a.) no longer 
factors. 


II. THE PROPAGATION OF ORDER IN TWO. 
DIMENSIONAL LATTICES 

As was already pointed out in the introduction, 
the characteristic value problem (1.4) for three 
dimensions, is very difficult to handle. Therefore, 
in the hope of obtaining useful qualitative 
information which might possibly be carried over 
into three dimensions, we turn to the two- 
dimensional case where the corresponding prob- 
lem is solvable. The matrix theory in two dimen- 
sions was developed mainly by Kramers and 
Wannier’ in a paper on the statistics of two- 
dimensional ferromagnets. 

In Section 3 we write down the matrix of 
Kramers and Wannier and review the properties 
of it which we shall find useful. Section 4 is 
devoted to a perturbation theory of the charac- 
teristic value problem, the main results of which 
are to show that the maximum characteristic 
value is doubly degenerate from T=0 to some 
finite temperature, and to obtain power series for 
the maximum characteristic value and the cor- 
responding characteristic vectors valid at low 
temperatures. In the remaining sections, the 
solution is applied to various problems of phys- 
ical interest. The long range order as a function 
of the temperature is obtained in Section 5, and 
in Section 6 a proof is given that an order- 
disorder transition actually occurs in the two- 
dimensional crystal. In Section 7 we discuss the 
general problem of the propagation of order 
through the lattice, showing in detail the gradual 
transition from the short range order to the 
influence at very great distances. In Section 8 
Zernike’s approximation is compared with the 
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Fic. 1. A two-dimensional crystal on the surface of a torus. 


exact solution. In Section 9, we treat the ordering 
influence of more complicated configurations of 
fixed atoms. 


3. The Matrix of Kramers and Wannier 


Kramers and Wannier have introduced a sim- 
plification into the theory of two-dimensional 
square lattices by making use of a special con- 
nectivity scheme for their crystal. This scheme 
can best be described as being equivalent to 
having the lattice sites regularly distributed 
along a continuous line twisting its way in screw- 
like fashion over the surface of a torus as in Fig. 1. 

m pitches of m atoms each in this case take 
the place of Z layers of MXN atoms each in the 
three-dimensional case. For convenience the 
atoms are distinguished by + and -- instead of 
A and B, coordinates 41, #mn Capable of 
assuming the values +1 and —1 being assigned 
to the respective lattice sites. The configuration 
of the atoms on a pitch will thus be given by a 
set of m values (n° - +241). which we denote by 
a for brevity. In the magnetic case » represents 
a spin capable of two orientations; for the case 
of alloys, the yu; are still convenient for describing 
the various configurations. The energies of 
interaction between neighboring pairs are written 
and V(+—)=}J, so 
that the total energy is expressible as 


ume, (3.1) 


where the sum is over all pairs of nearest neigh- 
bors. If J is positive, like atoms will attract one 
another and the ordered arrangement at low 
temperatures will be one in which the + and — 
atoms are separated into two phases. We shall 


OF ORDER IN CRYSTAL LATTICES 163 


then say that we are in the “ferromagnetic case.” 
The ‘‘antiferromagnetic case,’’ where J is 
negative, will lead to ordered states at low tem- 
peratures in which + and — atoms alternate on 
the lattice sites. 

As in the three-dimensional case, a matrix Jt 
is found whose characteristic value problem is 
fundamental for the statistical theory. Jt has 
2" rows and columns labeled by the 2” dif- 
ferent possible configurations of the atoms on a 
pitch, a, a’, a’, ---. However, the matrix is not 
symmetrical, and it thus becomes necessary to 
introduce right- and left-handed characteristic 
vectors, say A, and B,, where 


MA,=r,A, and B,M=r,B,. (3.2) 


With proper normalization the A, and B, can be 
made to satisfy orthogonality conditions. 


B,A, = By(a)Ag(a) = (3.3) 


It actually turns out that the analogue of the 
matrix (1.3) in the three-dimensional case, is not 
M, but the mth power of Mt. Thus one finds, upon 
expressing the matrix element (a:|9"|a2) as a 
bilinear combination of the components A,(a:) 
and B,(a;). 


2" a2) +0(a1, @ 


(3.4) 


where v(a2) is the total energy of interaction 
between the nearest neighboring atoms of pitch 2 
and v(a;, a) is the energy of interaction between 
the pitches 1 and 2. This relation, just as (1.5), 
plays a fundamental role. 

By repeated use of (3.4) and (3.3) one finds 
for the partition function, the expression 


(3.5) 


and it appears, by comparison with (1.7), that \" 
is the analogue of the \ of the three-dimensional 
case. The advantage of this result of Kramers 
and Wannier is that now \ has a meaning which 
is independent of the size of the crystal, the 
logarithm of Amax being essentially the free energy 
per particle. The variation of the properties of 
the crystal as it becomes infinite in two directions 
can therefore be studied conveniently by looking 
at Amax- 
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method chosen by Kramers and Wannier is ill I I 
the following: for any configuration, e.g.,  — | 3.7 
++——++-—+, one replaces every + by a0 v2 


and every — by a 1 (giving in our example 
00110010), reads the resulting number in the I being the identity matrix of order 2"~', we have 


dual number system and takes that as the order V.(K)| 0 . 
number of the configuration. For the configura- ||— 
tion above the order number is 50. It proves con- 7 0 | VK) » 38) 
venient to separate the configurations into two h ‘ 
classes in accordance with the sign of-za. The “°° 
configurations of the class u,=1 consist of those al 1 
with order numbers 0, 1, ---, 2*"'—1 and are 1 
arranged in that order. The remaining half of 4 
the configurations belong to the class = —1 and 
are arranged in the order 2"—1, 2"—2, ---, 2"-, V.(K)= a1 
The reason for the arrangement is the following. P 
Configurations a and & in corresponding places al . 
in the two classes have order numbers which add al 
up to 2"—1 and are conjugate to one another in 81. 
the sense that a can be obtained from & by 
changing every + to a — and vice versa. (3.9) 
Arranging the rows and columns to correspond eo 3 
to the arrangement of the configurations in the a 1 
two classes, and introducing the new variable K a 1 
and the functions a and 8 by 
1J V_(K)= 
=-—, B =e 
2 kT 
the matrix 2? has the form 
al | 
- | We shall call V,(K) the “plus V matrix” and 
V_(K) the “minus V matrix.’”’ The characteristic 
al | 00 values and characteristic vectors of It(K) fall 
00 | wl into two classes. If ¢,(K) and ¢_(K) are right- 
| handed characteristic vectors of V,(K) and 
| Bl V_(K), respectively, corresponding to the charac- 
| Bl teristic values \,(K) and \_(K), the associated 
(3-6) characteristic vectors of IM(K) are 
| al t,(K) t,(K) 
| A,(K)=H"™ - 
00 | al 0 t,(K) 
| 
| 0 q_(K) 
B1 | A-(K)=H" = 
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To write the matrix 22 as a square array, it is 
necessary to arrange the 2" configurations easily seen to be reducible. Thus, with a matrix 


a=(un-**uo1), in some definite order. The H defined by 


where blank spaces indicate zeros. It(K) is 
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(3.9) 


The components are arranged in the order of the 
two groups and wa=—1. Similarly, if 
$.(K) and ¢_(K) are the left-handed charac- 
teristic vectors of Vi(K) and V_(K), the cor- 
responding vectors for It(K) are 


B,(K)=||¢4(K), :(K)||, 


B_(K) =||¢_(K), 


As a final property of 90(K) which should be 
mentioned, there is a simple device by which one 
can treat the antiferromagnetic case once the 
ferromagnetic case has been discussed. In the 
former case, the interaction energy J is negative 
and the parameter K is therefore negative. In 
the ferromagnetic case, K >0. If R is a permuta- 
tion matrix which changes every alternate atom 
on a pitch from + to — and vice versa, it is 
readily verified that 


RBW(—K)H“R™ 
Vi(K) | 0 


0 | —V_(K) 


so that the characteristic values of J(—K) are 
\,(K) and —\_(K). The characteristic vectors 
A,(—K), etc., are obtained from A,(K), B,(K), 
-++ by permuting the components according toR. 

The problem with which we are now faced is 
that of finding the characteristic values and 
characteristic vectors of the plus and minus V 
matrices as functions of the temperature. 
Kramers and Wannier have neglected the minus 
V matrix in their treatment, but this is unim- 
portant for the discussion of the thermal proper- 
ties. We are, however, interested in the de- 
generacy of \max and shall find in the next section 
that the maximum characteristic value of V_ is 
equal to the maximum characteristic value of V,. 
in the limiting case of an infinitely large crystal. 


(3.11) 


» (3.12) 


4. Solution of the Matrix Problem for 
Low Temperatures 


For the low temperature ferromagnetic case, 
K is very large and positive, so that the param- 
eter 8=e~** is very small, approaching 0 as the 
temperature tends to zero. It then becomes 
possible to carry out an expansion in powers of 8. 
Thus, by defining a new matrix, U,(8), the plus 
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V matrix can be written 


18 
18 
U..(8) can be expanded in the form 
10 
10 
U,(8)= 
01 
01 
‘ol 
01 
01 
o1 
01 
+6 (4.2) 
10 
10 
10 


Defining the constant matrices Uo, U:, and U, 
by (4.2) we may write for brevity 


U,.(8) =Uo+601+ 6°02. (4.3) 


If there are m atoms per pitch, these matrices will 
have 2"—' rows and columns corresponding to 
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the configurations 0, 1, 2, ---, 2"-'—1 which all 
have 

Our procedure will consist in obtaining power 
series in the parameter 8 for the maximum 
characteristic value Amax and corresponding 
characteristic vector Umax of the matrix U,(). 
It is evident that U> has a single non-degenerate 
characteristic value of unity and a 2"-'—1-fold 
degenerate value of zero. Therefore, Amax will be 
well separated from the other characteristic 
values at low temperatures and will be an ana- 
lytic function of 8 in the neighborhood of B=0. 
One will thus be justified in trying series ex- 
pansions for Amax and Umax of the form 


Y= tax = 


Inserting (4.4) into the equation U,(8)¢=Ad, 
and equating coefficients of like powers of 8, we 
find a set of equations for the \, and ¢,: 


0) Voto =Aoho, 

1) (4.5) 


(4.4) 


The characteristic vector t%» corresponding to 
Xo = 1 evidently has only one non-vanishing com- 
ponent, that corresponding to configuration 
number zero. We therefore take Yo(0)=1 and 
tho(a)=0 for a=1. 2, ---, 2" "*—1. Thus for vo 
and its adjoint vector, to, we have 


0 

Y=]10--- Of. (4.6) 

0 

Since U:to=0, the equation for and is 
= +Ardo. 


Taking the inner product on the left with the 
vector to and noticing that YoUo= tho, we find 
that }4,=0. The vector %, then satisfies the 
equation Uot:= q:. It is clear, however, that we 
are at liberty to choose q; and all the remaining 
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t, orthogonal to to, saving the final normaliza- 
tion of ¢ until the end. With this restriction, 
1 must vanish identically. 

_ Multiplying the equation for q2 on the left by 
to, we find in this case also that A2=0. To find 
2 we must therefore solve the equation 


(I—Uo) = 


where I is the identity matrix. 

Denoting the components of 2: by (a) 
(a=0, 1, 2, ---, 2"-'—1), we find for (I—U,)y, 
and Us the vectors 


4(1) 
4 (2) — 
4(j) — 4(23) 


— (2-1-2) 


0 
0 


0 
0 


0 
1 


Equation (4.7) is of the type (I—U»)¢=, 
where & is a known vector and ¢ is to be found. 
The entire perturbation procedure consists in 
repeatedly solving equations of this kind. The 
method used is the following. Since the com- 
ponents ---, of (I-Uy)¢ 
are identical with those of 4, they are obtained 
immediately from the corresponding components 
of &. The remaining half of the components of ¢ 
are obtained with the help of those already 
known. For example, 


4(2"-*—1) = + 1), 
2) = —4) + 
etc. 
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A component of even index j is coupled di- 
rectly to the component of index $j and indi- 
rectly to those with indices }j, $j . . . if jis also 
divisible by 4, 8. . . . One continues to follow 
the even components j/2* until further division 
by 2 is no longer possible. 

Applied to t2, the procedure om — —1) 
=1 and for a=0, 1, 2, » 
Multiplying the equation for ¥; on ‘ie left by 
to, we find \s= toUite. In general, from the fact 
that % is orthogonal to the other ¥; and that 
oU:=0, we find Ar = However, since 


%oUi=|/0 1 0 --- Ol, 


we have 


t,-a(1). (4.9) 


Since the component number 1 of g. vanishes 
we see that As=0. ws is then to be found from the 
equation (I—Uo)¥s= The non-zero com- 
ponents of Uitf2 are 


(Uitte) (2"-*—1)=1, (Cite) (2"-*) =1. 
The non-zero components of 3 are therefore 


q3(2")=1, =1, 
t3(2*) =1 (k=1, 2, +++, n—3) 


arising from the connection with the component 
2*-*, and ¥3(1)=1 arising from the connection 
with §3(2). 

From (4.9) we then find 4=3(1)=1. The 
vector t, therefore satisfies the equation 


(I— Uo) = Vids — to. 
However, Since 
(Wig) (7) = (Tig) 


~1-§)=4(2j+1), 


we find for the non-vanishing components of 
the following: 


(Wis) (0) = (Wigs) — 


1)=1 


(Wigs) (2"-* — 1) = (Wigs) (3.2"-*) = 1. 
Therefore, 
— 
and from this 
(k=1-- 


finally 44(2"-*—1) =1. 


1)=1, 


-m—4), t4(3)=1; 
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This perturbation procedure was carried as far 
as the 12th power of 8 and gives for the maximum 
characteristic value of the plus V matrix" 


1 
Ay + 26° +56%+ 148" +4482 + 4.10) 


The convergence of this series will be inferred 
later on. More complete series for 4, will also be 
given later together with the characteristic 
vectors of the minus V matrix. However, up to 
terms of the order 8°, the non-zero components 
of the approximating vectors are given in Table I. 

The expression (4.10) for A,, if carried out 
completely, would essentially give the free energy 
per particle for an infinite crystal (i.e., one for 
which m=). For any finite crystal the cor- 
responding series will agree with the limiting one, 
(4.10), up to a certain power of 8, after which 
they will differ. It is not difficult to see from the 
perturbation procedure that the agreement will 
persist to higher and higher powers of 8 as n 
increases. The only reason for a difference 
between the limiting series and the one for 
finite m, is that at some stage in the calculation 
two or more components of a vector which are 
labeled differently actually become identical. For 
example, in the case m= 5 the component number 
3 of the fourth approximating vector q, does not 
have the value 1 as one would find by the above 
methods, but actually has the value 2. This 
comes about because the component labeled 
2*-*—1 is also the component number 3. Because 
of these coincidences the components of the 
approximating vector §, and hence also the coef- 
ficients 4, for a given n, will eventually differ, as 
r increases, from the &, and 2, for larger values 
of n. However, these coincidences can be post- 
poned by taking larger and larger values for n, 
so that at any fixed stage r in the approximation 
a given component of ¥, or the value of 4, will 
be the same for all sufficiently large finite 
problems. 

The same type of perturbation procedure is 
applicable to the problem of the minus V matrix. 
Carried as far as terms of the order 6" we find 
that its maximum characteristic value \_ is 


“4 The terms up to §"° of this series were | 
tained by Bloch, method. od. FB 
f. Physik 61, 206 (19. 
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TaBLe I. The approximating vectors th to ts. 


Value} 


1 


2* k=1---n— 

7.2% k=1--- 

20-1 


identical with \, for sufficiently large . More- 
over we find to this’ approximation that the 
components of the vector t_ are identical with 
those of 4 except for characteristic differences 
in sign. For any finite value of m, the charac- 
teristic values \, and \_ will differ. However, 
experience with the perturbation method leads 
us to believe that as m becomes increasingly 
large, the difference between A, and \_ is post- 
poned to higher and higher powers of 8, so that 
in the limiting case of an infinite crystal the 
degeneracy between \, and \_ will be expected 
to hold as far as the series converge. Really to 
give a proof of the degeneracy, one would have 
to devise some inductive argument of the type 
given by Lassettre and Howe'® to show the 
degeneracy of the maximum characteristic value 
of the matrix for the three-dimensional problem. 
_ To find the left-handed characteristic vectors 
B, and B_ the problems of obtaining the left- 
handed characteristic vectors of the plus and 
minus V matrices are treated by the same per- 
turbation methods. The same series for 44 and A_ 
are found as previously and the equality, except 
for characteristic sign differences, of the com- 
ponents of ¢, and ¢_ is also observed. 

For a proof of the convergence of the series for 
A, and A_ we refer to the work of van der 
Waerden.® By a consideration of the probabilities 
for finding polygons of various lengths enclosing 
atoms which are not ordered with respect to one 
another, he calculates the long range order and 
the energy of the crystal at low temperatures, 
obtaining series identical with ours. Although 
his method would be very cumbersome for 


4% E. N. Lassettre and J. P. Howe, J. Chem. Phys. 9, 
801 (1941). 


treating the general propagation problem, it is 
more powerful than the matrix method in the 
respect that it provides simple convergence 
proofs. From the convergence of the series for 
the energy we can conclude that the series for 
d, and \_ also converge. The energy is obtained 
from the partition function, or from the partition 
function per particle, by 
8 log Z 7 d log Ay 


0(1/kT) dg 
d log (8+) 


=—NJ+NJ 


= —NJ+NJ{48'+ 126°+ 368" 
+1208" ---}. 


This is the same as the series for the energy which 
is given by van der Waerden and from his proof 
of its convergence we see that the logarithmic 
derivative of BA, is analytic in the neighborhood 
of 8=0. The logarithm of 8d, is therefore ana- 
lytic at 8=0 and hence the same is true of BA, 


itself. 
5. The Long Range Order 
The probability P(a:,a,) of finding simul- 
taneously pitch 1 in configuration a, and pitchr 
in configuration a, is given by an expression 
analagous to (2.1) of Section 2. Using (3.3), 
(3.4), and (3.5) one finds 


P(a, a,) = La Ap 


mn 
P 


x B,(a:)A,(a 1)B,(a,)Ap(ar) 


In the ferromagnetic case, for sufficiently low 
temperatures, we have found the maximum 


(5.1) 
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characteristic value to be two-fold degenerate. 
Therefore, for two pitches far enough apart so 
that r—1 and m—r-+1 are large integers, we 
have 


P(a1, ar) = 31 (ar) +B_(a1)A_(a,) } 
X (a,)A(a1) (5.2) 


Let us now ask for the probability of finding 
two + atoms simultaneously, one in a given 
position in pitch 1 and the other at a given place 
in pitch 7. In particular we might ask that the 
two end atoms of the separate pitches be both +. 
Then the probability P in question is obtained 
by summing P(a:, a,) over those configurations 
consistent with having + atoms at the end of 
the pitch. These comprise the first half of all the 
configurations, namely, the ones numbered from 


Oto 21-1. 


From (3.10) and (3.11) the normalized forms 
of A, and B, are 


Ay 
(5.3) 


By = (29445) 


In the same way, the normalizing factor for the 
minus vectors is (2¢_4_)—4. However, the series 
expansions show that ¢,4,. It also turns 
out that B,(@)A;(@)=B_(a)A-(a) and that 
B,(a)A_(a) = B_(a@)A,(a@). Using these relations, 
multiplying out (5.2), and summing over the 
first half of the configurations, we obtain for P 
the expression 


(S.4) 


It can be verified without very much difficulty 
that this same expression remains valid if the 
two + atoms under consideration do not neces- 
sarily occupy corresponding positions in the two 
distant pitches. P is therefore the chance of 
finding simultaneously two + atoms very far 
apart in a ferromagnetic crystal. In the absence 
of long range order, P= }. 

If the crystal is antiferromagnetic, this prob- 
ability will of course depend on the odd or even 
character of the number of steps between the 
two atoms. Instead of a degeneracy in Amax, We 
have, from (3.12), 4,(—K)=—A_(—K). The 


total number of atoms in the crystal, mn is even, 
however, since there are as many + atoms as 
— atoms. Thus, \,(—K)""=\_(—K)™*, and the 
effect is the same as though there were a real 
degeneracy. Also, from (3.12), we see that the 
characteristic vectors ¢1(—K), are 
obtained from ¢+(K), ¢4(K) by permuting the 
components in accordance with a certain per- 
mutation matrix R. The components of ¢+(K) 
and %+(K) are permuted in the same way, 
however. Consequently the inner products of 
the vectors are invariant; that is, for example, 


$.(K)-q-(K) = 
= 


Carrying out the computation for P, we find an 
expression 


(5.5) 


where & is, say, the smallest number of steps 
between the two + atoms entering into the 
definition of P. If k is even, the two + atoms form 
an ordered pair, and P exceeds }. If k is odd, the 
pair of + atoms is not ordered, and P is less 
than }. 

To make connection with the usual long range 
order S we shall obtain an expression for P in 
terms of S and compare with (5.4) and (5.5). 
For an antiferromagnetic alloy, one divides the 
sites of the crystal into two groups, a and 8 sites, 
the former being occupied exclusively by + 
atoms at the zero of temperature and the latter 
by — atoms. At a finite temperature one denotes 
by r the probability that an atom be right, or on 
a site which is appropriate to it, and by w the 
probability that a + or a — atom be wrong, 
that is, be on a site not appropriate to it. The 
long range order S is then defined by writing w 
and r in the form 


w=$(1—S), r=3(1+5). (5.6) 


In terms of S, the probability, P, for finding 
two atoms, of which one is +, and which are 
correctly ordered with respect to one another is 
given by 


+ —) (5.7) 


and 
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This can be seen as follows. If one takes the 
factors in order, the first } is the chance that an 
a site is selected for the first atom of the pair, 
which we may take as +. The chance of getting 
a + atom on this site and also of getting the 
right atom at the distant site, is then }(1+.S)’. 
If the site chosen originally is a 8 site, the two 
atoms of the pair would both have to be wrong 
from the point of view of the a—8 sites in order 
to be correct with respect to one another. This 
would give the factor }(1—.S)?. In the same way 
the probability for finding two distant atoms 
which are not correctly ordered with respect to 
one another, is }(1—.S*). 

From (5.4) and (5.5) it therefore follows that 
the long range order S is given in terms of the 
characteristic vectors by 


S= 


The computation of S up to terms of order 8° 
is given in Table II. All configurations a which 
give the same series expansion for ¢4(a)¥_(a), 
to terms of order 8*, are grouped together and 
listed in the first column. The number of con- 
figurations in each group is given in column two. 
The third column gives the coefficients of 6°, B*, 
8°, and 6* appearing in the series for ¢4(a)t_(a). 
To obtain ¢,4_, we must sum ¢,(a)§_(a) over 
all a from 0 to 2"-!—1. This is conveniently done 
by multiplying the coefficients in the third 
column by the number of configurations in the 
second column and then adding. The results are 
indicated in the last column. We have 


(3n—7) 


(5.8) 


1 17 
2 2 


The product ¢,%, is gotten in the same way. 
Its table differs from that of ¢,4_ only in that 
all the signs of the ¢,(a)d,(a) are positive. 
Otherwise the numerical values are the same as 
for the ¢4(a)$_(a). Making the appropriate 
changes we find 


$44, =1+n8'+ 


1 21 
+ (6.90) 
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Both ¢,4_ and ¢,4, depend on the size of the 
crystal, m, but their quotient, which is the long 
range order, must be independent of m for very 
large n. Performing the division, we obtain to 


terms of order 6” 
= 
(5.10) 


and 
P=}(1+S") =} — 168° 
— 688" 


(5.11) 


This same series for S was found by van der 
Waerden and a proof was given that it is con. 
vergent for 8<4. The question of the actual 
radius of convergence of the series is still open, 
however. A highly probable conjecture can be 
made in the light of the investigations of Kramers 
and Wannier. On the basis of an important sym- 
metry property of the plus V matrix, they are 
able to show that if it is assumed that the crystal 
can have only one point of singular behavior, 
that point must be given by 


B.=Vv2—1=0.414. 


(5.12) 


As we shall prove in the next section, there is no 
long range order at high temperatures. It there- 
fore seems very probable that the point 8,= 0.414 
represents the Curie point for the transition from 
order to disorder and that the series for S has 
radius of convergence 0.414. A real proof of this 
fact and an investigation of the limiting value of 
S(8) as 6 approaches 8, from below, would be 
very much desired, however. 


6. Proof That an Order-Disorder Transition 
Actually Takes Place 


Although the location of the transition point 
between states with finite long range order and 
those with no long range order has not been 
rigorously obtained by the matrix method, it is 
easy to give a proof that such a transition really 
occurs. Such proofs have been given by Peierls" 
and van der Waerden® using the properties of 
polygons enclosing regions of disordered atoms. 


1 R. Peierls, Proc. Camb. Phil. Soc. 32, 477-481 (1936). 
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TABLE II. The calculation of ¢,1-. 


$+ (a) 
Configuration, : 

1 1 
k=1---n—2 n—2 § n—2 2n—4 5n—10 
3.2" k=1---n—3 n—3 1 6 n—3 6n—18 
2i+2* ju3--+-n—2, R=1-++j—2| }(n—4)(n—3) 1 jn*— 6 
7.2* k=1---n—4 n—4 1 n—4 
1 —4 
1 1 1 8 
3 1 . 1 7 
7 1 1 1 
2k+1 k=2---n—2 n—3 1 n— 3 
1 -1 
1 -7 -1 -7 
n —n+3 
1 -8 -1 —8 

n—2 


From the point of view of the matrix theory all V_(0) have the forms 
that is necessary is to show that for sufficiently 


high temperatures the maximum characteristic ee rs 
value of I2(K) is no longer degenerate. 
For high temperatures K(=}J/kT) tends to 
zero and both a@ and tend to unity. Thus for 
very high temperatures the matrices V,(K) and 
111 --- 111 
V_(K) are practically 111 


11 (6.2) 
11 vas 


V..(0)" consists entirely of ones while V_(0)’ has 
‘4 alternate rows of ones and minus ones. Since the 

11 matrices are of rank one, their characteristic 
equations reduce to the form 


V_(0)= 1 1} (trace (6.3) 


sai The trace of V,(0)’ is 2”. Therefore the charac- 
—-1i-1 teristic values of V(0)’ are 2’, 0, 0, ---, 0 and 
hence the maximum characteristic value of V(0) 

is 2 and all the rest are 0. The trace of V_(0)’ 
Taking a finite problem for which the matrices vanishes, however. As a result all of its charac- 
are of the order 2", one finds upon successive teristic values and hence all the characteristic 
multiplication that the rth powers of V(0) and values of V_(0) vanish. 
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This proves that for sufficiently high tem- 
peratures the maximum characteristic value of 
$M is non-degenerate and that there is therefore 
no long range order. 


7. The Propagation of Order 


The principal advantage that can be claimed 
for the matrix method is that it provides a solu- 
tion not only to the problem of the long range 
order, but also to the more general problem of the 
propagation of order throughout the lattice. For 
example, suppose that a certain site in the 
crystal is known to be occupied by a + atom. 
We can then investigate how this knowledge 
affects the probability for + (or —) atoms at all 
other sites and see the gradual transition from 
the strong short range influence of the fixed + 
atom to the weaker long range influence which 
is described by the long range order. This is the 
problem originally proposed by Zernike. More- 
over, using the matrix method we can investigate 
the propagation of order for more complicated 
initial configurations consisting of several atoms. 
Thus we shall find the ordering effect of a pair 
of fixed atoms, of combinations of three, four 
atoms, and so on, in each case observing the 
change from the short to the long range effect. 


Zernike’s Problem. The Transition from Short to 
Long Range Order 


If the expression (5.1) for P(a:, a) is summed 
over all configurations a,, one finds for the 
probability of a configuration a, in pitch 
number 1 


P(a)= 
Dip 
At low temperatures and for ferromagnetic 
crystals, the twofold degeneracy of Amax leads 
to the expression (dropping the subscript 1) 


P(a) = 3(B,(a)A, (a) +B_(a)A-(a)) 
=B,(a)A,(a). (7.1) 


This is in marked analogy with the quantum 
mechanics, where the probability density is pro- 
portional to the wave function multiplied by its 
complex conjugate. From (7.1) we can also see 
that the theory is entirely symmetrical in the 
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two kinds of atoms, + and —. Thus if we make 
a transformation which takes a configuration q 
into one & obtained from it by changing all + 
atoms to — atoms, and conversely, we find 
from (3.10) and (3.11) that B,(a@)=B,(@) and 
A,(a)=A,(@), so that P(a)=P(a). This sym. 
metry is evidently necessary, since + and ~— 
atoms have been treated equally from the start, 
the energy of interaction between a pair of + 
atoms being equal to the energy of interaction 
of a pair of — atoms. 

Any question concerning the probability of a 
certain disposition of atoms on lattice sites all in 
one pitch is answered by summing P(a) over all 
configuration a@ consistent with the restrictions 
proposed. We shall, therefore, always have to 
deal with the expressions of the type 


$+(a) t,(a) ™ -V 


264%, 


By (a)A,(a) = 


where the prime indicates that the summation 
is to be taken over the configurations allowed by 
the restrictions, and V is the sum of ¢+(a)¥+(a) 
over the suppressed configurations. The sum V is 
generally easier to carry out than the direct sum 
over the allowed configurations. 

Let us now apply (7.2) to the investigation of 
the ordering influence of one atom on its neigh- 
bors in the same line of the crystal. Let P(+ 1), 
«++, P(+uir,)," be the probability of 
finding a + atom at a site conveniently called 1 
and then an atom rightly ordered with respect 
to it at a distance of one site away, two sites 
away,... , (p—1) sites away, respectively. We 
can write P(+.7,) as equal to the probability 
for finding a + atom at the site 1 in the lattice 
multiplied by the probability that the right atom 
is found p—1 sites away when it is already known 
that a + atom exists at site 1. The latter corre- 
lation probability we denote by C(+:|r,) and 
take it as a measure of the ordering influence of 
a single + atom. Of course 


P(+1ry) 


17 This r notation is used so as not to distinguish between 
the ferro etic and antiferromagnetic cases. The calcu 
lations will be made for the ferromagnetic case, but it will 
be shown that the results are valid in the other as well. 
For the ferromagnetic case, r2, for example, would be +s, 
while for the antiferromagnetic case it would be —:. 


C 
C 
an 
mi 
tet 
(8 
fre 
sit 
1 the 
inc 
the 
lim 
enc 
the 
tha 
F 
C(- 
cor 


These correlation probabilities are easily cal- 
culated from a table of values of ¢4(a)t,(a) 
similar to that given in Table II. For example, if, 
in (7.2), we make the restriction that the last two 
atoms, numbers m—1 and m—2 are both + we 
have the conditions necessary for the calculation 
of C(+:|r2). Similarly, for C(+1|rs) we would 
hold the atoms m— 1 and n—3 fixed as +. We find 


C(+:|r2) = 1 — 66° — 188° 
— 606" — 
C(+1| 17s) = 1—26* — 86° — 3088 
C(+1| 174) = 1 — 28* — 86° — 3288 
nin 
4 
C(+:| 15) = 1 — 26*— 86° — 326° 


C(+1|r) =1—26'—86*— 326° 
— 1368" —6148"—--., 


C(+1| rx) = 1 3288 
— 1368" — 61482 


and C(+:|r,)=C(+:|rs) to terms of order 8” 
for p=7. 

The correlation probability, C(+.:|r,), has its 
maximum value of unity at the absolute zero of 
temperature (8=0) and for finite temperatures 
(8>0) falls off monotonically as one goes further 
from the central site at which the + atom is 
situated. Examining the sequence (7.4) we find 
that the coefficients of the various powers of 8 
increase monotonically in absolute value until 
they reach a certain limit, the approach to the 
limit being rapid at first and then slower at the 
end. For the lower powers of 8 this approach to 
the limiting value occurs earlier in the sequence 
than it does for the higher powers. Thus, from 
(7.4), 


C(+1{1r2) > C(+1| 17s) > C(+1| 14) 
(7.5) 


For p=6 we see that to terms of the order 6", 
C(+:|r,) becomes identical with the long range 
correlation (5.11). Denoting the long range cor- 
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relation by C(+:|r.) we have 
C(+:1| ro) = +S). (7.6) 


In this approximation the long range effect is 
already attained five sites away from the central 
atom. This holds of course only for those tem- 
peratures for which terms of order higher than 
B can be neglected. As the temperature in- 
creases the transition to the long range limit 
takes place only after one goes considerably 
further from the central site. 

The correlation probability C(+,|r2) is closely 
related to the average energy of the crystal. Thus 
if Q, stands for the probable number of right 
pairs of atoms each having an energy of inter- 
action —4$J and Q, for the number of wrong 
pairs of atoms each having an energy of $/, we 
have for the average energy 


E=- 


Remembering that for a crystal containing NV 


atoms there are 2N pairs of nearest neighbors, 
we find x 


and 

= —NJ+NJ(48*+ 126°+ 366% 
(7.7) 


This series for E has been given before in (4.11) 
and was obtained there from the expansion of 
the maximum characteristic value “A,. The 
present method is based on the characteristic 
vectors and the agreement between the two 
provides a check on the internal consistency of 
the results. 

It is not necessary to confine oneself entirely 
to atoms on the same line with the fixed + atom. 
For in virtue of the screwlike construction of the 
lattice, the atom at the end of the pitch, »—1, 
is situated near the atoms 0, 1, 2, --- at the 
beginning of the pitch, as shown in Fig. 2. We 
can find the influence of a fixed + atom at sites 
on the next line by applying (7.2) with the re- 
strictions that atoms m—1 and 0, m—1 and 1, 
n—1 and 2, etc. are to be successively held 
fixed as + atoms. The notation C(+'|r,) suggests 
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Fic. 2. The m sites of a single pitch. 


itself for the correlation.probabilities. The results 
are 


C(+*|r2) = 1 — 284 — 88° — 2688 
— 888" — 3266" 

1248" —4g4g"—..., 

= 1 —264— 8° — 3288 
~1366"— 

C(**| = 1 — 284 — 88° — 3268 
1368" 

C(**| = 1 —26*— — 32,8 

1368" — 6148"... 


(7.8) 


This sequence has the same properties as the 
sequence (7.4). A comparison of the two shows 
that 


C(+1|1r2) > C(**|1r2) > C(+1| 13) 
> C(+*| > C( +1] 174) > 14) > (7.9) 


indicating again a monotonic decrease of the cor- 
relation probability. as the distance from the 
fixed atom increases. 

The correlation probabilities (7.4) and (7.8) do 
not constitute a complete solution of Zernike’s 
problem since there are many distances between 
pairs of atoms which are not included. However, 
this is as far as one can go without a knowledge 
of the entire spectrum of characteristic values of 
the V, matrix. This is because the expression for 
P(a, a), where r is a small integer, cannot be 
simplified to the point where only the maximum 
characteristic value of V, enters. Nevertheless 
the range of distances covered by (7.4) and (7.8) 
is fairly complete and for those that are left out, 
one can get good estimates of the coefficients of 
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the smaller powers of 8 by applying the condition 
for monotonic decrease with distance. 

It should be pointed out that the expressions 
found for the correlation probabilities are valid 
for antiferromagnetic crystals as well. As we see 
from (3.12), the vectors Bs(—K), A+(—K) 
applicable in this case have components which 
are identical except for the order with those of 
Bs(K) and As(K) appropriate to the ferromag- 
netic case. The permutation matrix R, carrying 
the one set of vectors into the other, carries each 
configuration into the one obtained from it by 
changing the sign of every alternate atom in it, 
Thus for C(+:]| +2), calculated in the ferromag- 
netic case, we summed over all the configurations 
for which atoms »—1 and n—2 were both +, 
The permutation matrix R carries these con- 
figurations into all those for which the atoms 
n—1 and m—2 are, respectively, + and -, 
Therefore C(+:1]++2), as calculated in the ferro- 
magnetic case, is equal to C(+,|—2) for the 
antiferromagnetic case. 


8. Comparison with the Results of Zernike 


Carrying out Zernike’s method for treating the 
propagation problem, we find in the case of two 
dimensions for the long range order and the cor- 
relation probabilities 


S=1—26*—8* — 2088 
C(+1| 72) 
C(+| 72) = 1 
C(+1|rs) 
C(+1| rs) 
PES (to terms 
of order 


(8.1) 


Comparison with (5.10), (7.4), and (7.8) shows 
that Zernike’s approximation at low tempera- 
tures is valid up to terms of the order 6°. 

In three dimensions we make use of the series 
for S and C(+:|r2) found by van der Waerden.! 
The comparison with Zernike is given in (8.2). 


Vi 


174 
a-3 n-2 ° 

Z 
Ze 
in 
er. 
di 
th 
ne 
pe 
ate 
inc 
acl 
nit 
col 
| inf 
din 
cer 
We 
get 
din 
lea: 
| is 
rivi 
9." 
N 
arré 
adje 
prol 
and 
at si 


A) 


van der Waerden: 
S=1—26*—128"-+ 148" —908" 


Zernike : 
van der Waerden: 
C(+:|r2) =1—26*— 148" 


— 708" + 1768'S — 


(8.2) 


Zernike : 
C(+1| 72) =1—26°+668+ ---. 


Zernike’s method therefore seems to be better 
in two dimensions than it is in three. This can 
perhaps be understood from the following consid- 
eration. The essential approximation in Zernike’s 
difference equation comes from the assumption 
that the probabilities for the atoms forming the 
nearest neighbors of a given atom, are inde- 
pendent of one another. Although no two of the 
atoms at these surrounding sites are nearest 
neighbors, the probabilities for them are not 
independent because of the possibility of inter- 
action through intermediate atoms. The mag- 
nitude of the correlation will in some way be 
connected with the number of short paths of 
influence between the two atoms. In three 
dimensions, the correlation will presumably be 
greater because the number of available paths is 
certainly much greater than in two dimensions. 
We might therefore expect the approximation to 
get poorer as the dimension is increased. In one 
dimension, where the assumption of Zernike is 
least serious, it can be shown! that his solution 
is actually exact, although the method of de- 
riving the equation involves an approximation. 


9. The Ordering Influence of More Complicated 
Configurations 


Next to a single + atom, the simplest initial 
arrangement to consider is an ordered pair of 
adjacent atoms, +72. Denoting by P(+rer;) the 
probability for finding a plus atom at a given site 
and two atoms rightly ordered with respect to it 
at sites adjacent and (j—1) steps away, we have 


#8 A. Nordsieck (unpublished). 
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for the correlation probability 
C(+are| 13) = (+12). 
Application of the methods of Section 7 gives 


C(+1r2| 1s) = 1—6*—46*— 1768 


C(+1re| = 1 58% — 2488 
— 1098" 


C(+1r2| rs) = 258° 


C(+1re| 76) = 1 25,8 

— 1228" — 599g" ..., 
C(+are| rz) = 1—6*—Sp* — 

— 


and to terms of order 8” 
p=s. 


Comparison of (7.10) and (7.4) shows that in 
general 


3) >C(+1| 75-1). 


Thus, as expected, an ordered pair of atoms is a 
more effective ordering influence than a single 
atom. The sequence (9.1) shows the same kind 
of approach to the long distance limit as do the 
sequences (7.4) and (7.8). 

It will be interesting to compare the correlation 
probabilities for the ordered pair +72 with those 
for the ordered pair +'r, where the atoms are 
separated by a distance v2 times the lattice 
distance. (See Fig. 3.) By imposing, in (7.2), the 
condition that the atoms »—1, 0, and j—2, of a 
single pitch be held fixed as + we find the prob- 
ability P(*rer;). Dividing this by P(+r2) gives 


e 
+ 
3 4 j 


Fic. 3. An ordered pair of atoms separated by a distance V2 
(in units of the lattice distance). 
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Fic. 4. The shortest path of influence from site 1 to site 5 
when site 2 is occupied by a known atom. 


the correlation probability C(+r2|r;). The result 
1S 

173) =1—B*— 36°— 148° 
— 628" 

C(*r2| 74) = 1 —B*— 48° — 2088 

C(+r2| 75) = 218° 

C(*#r2| 76) = 216° 

C(*r2|1r7) =1—p*— 46° — 2188 


(9.2) 


and to terms of order 8" 


The comparison of (9.1) with (9.2) shows that 
the pair +72 is a more effective ordering influence 


than the pair +172; that is 
C(#72| 175) > 


The reason for this can be easily understood. In 
the diagram of Fig. 4 let us consider the corre- 
lation probability at site 5, perhaps, for an 
ordered pair of atoms at sites 1 and 2 and then 
for an ordered pair at 1’ and 2. With the atom 
at 2 held fixed, the shortest path of influence 
from 1 to 5 consists of six steps and is shown by 
the dotted line. However, for the atom at site 1’, 
the shortest path of influence consists of only 
five steps. Thus, if site 2 is occupied by a known 
atom, placing a properly ordered atom at site 1’ 
will be more effective at 5 than placing an ordered 
atom at site 1. We should, therefore, expect 
C(+'r2|1r5) to be larger than C(+72|rs). 

If the correlation probabilities arising from an 
ordered pair of atoms are greater than the cor- 
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relation probabilities arising from a single atom, 
we should find that the correlations for a dis. 
ordered pair will be smaller than those for a 
single atom. We shall denote such a disordered 
pair, or “dipole,” by +1w2. The correlation 
probabilities C(+1we|r;) can be obtained from 
results already given. Thus, expressing the 
identity 


P(+ wor j) +P(+srer = 
in terms of correlation probabilities, we have 
C(+1| 173) — C( +1] 172) C( +172! 75) 
1—C(+:|r2) 


= 
(9.3) 


Here r; means right with respect to the + atom 


at site 1. 
From (7.4), (9.1), and (9.3) we have then 


 — 
C(-+ 7s) =4(1 — 


(9.4) 
C(+1w2| 76) = 3 (1 — f—---), 
C(+1w2| =F (1 


and to terms of order f° 
p27. 


These series show that at large distances the 
correlation probability arising from a dipole, 
+ w2, approaches the value 4. The dipole there- 
fore has no ordering influence at all at large 
distances, since the probability of finding a + 
or a — atom at a given lattice point when the 
dipole is not present is also equal to 4. The atoms 
of the disordered pair thus have the tendency to 
annul each other’s ordering effect, the cancel- 
lation being complete at infinite distance. 

Many further examples of sequences of cor- 
relation probabilities could be given. Instead, 
however, we shall give in Table III a list of 
correlation probabilities for a line of the lattice 
from which all other correlations can be obtained 
by suitable combination. The table gives the cor- 
relation probabilities as a function of position in 
a line of the lattice for a sequence of initial con- 
figurations beginning first with one fixed + 


.3) 


4) 


F+9 F 


28.5 


+8588 


atom, then an ordered pair of atoms, then 3, 4, 
5, 6, etc. . . . consecutive ordered atoms. 

For each fixed initial configuration, the cor- 
relation probability decreases monotonically with 


TABLE III. Correlation probabilities for the atoms in one 
line of the lattice. 


1 -—-2 -6 -—-18 — 6€ —224 
C(+:\ 73) 1 —30 -—112 -—430 
C(+:1 1 —-2 —32 -—134 —576 
C(+:| rs) 1 —8 -—136 —612 
C(+:1| 16) 1 -—32 -—136 -—614 
C(+:|12) 1 —-2 —32 -—136 —614 
C(+1r2| 7s) 1 —-1 —4 -—-17 — 70 —291 
ra) 1 —-1 -—24 -—109 —487 
rs) 1 -—-1 -—25 —121 —580 
C(+1r2| fe) 1 —-1 —5 —25 -—122 -—599 
C(+1ra| 17) 1 —-1 —25 —122 —600 
1 —-1 —122 —600 
1 —-1 — 50 —241 
rs) 1-1 —4 -17 — 78 —376 
1 —-1 —4 — 89 —451 
12) -18 — 90 —469 
rs) 1 —-1 -—-18 — 90 —470 
C(+ rs) 1 —4 -—18 — 90 —470 
C(+1- + rs) 1 -3 -10 — 38 —170 
1 -16 — 65 —286 
C(+1- + 77) 1-1 —4 —-17 — 76 —360 
C(+1° 78) — 77 —378 
C(+1-+ 79) —4 -17 — 77 —379 
C(+:1-+ +14! rio) 1-1 —4 —-17 — 77 —379 
76) 1 -1 -10 — 37 —150 
77) 1-1 -16 — —265 
C(+1-+ 1-1 —4 -17 — 75 —339 
C(+1-+ 79) 1-1 —4 —17 — 76 —357 
C(+1° + 711) 1 -1 —4 -17 — 76 —358 
C(+1°+ 1 -1 -—3 -10 — 37 —149 
C(+1- + 78) 1-1 —4 -16 — 64 —264 
179) 1 —-1 —4 —-17 — 75 —338 
C(+1- + 710) -17 — 76 —356 
+ 711) 1-1 —4 -17 - 76 —357 
C(+1- 712) 1-1 -—-17 76 —357 
C(+1° + 18) 1 -—3 -—-10 — 37 —149 
19) -—4 -16 — 64 —264 
P10) 1 -1 —4 -17 — 75 —338 
| 1 —4 -17 76 —356 
112) 1 —4 —17 — 76 —357 
—4 -17 — 76 —357 
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distance and tends to a limiting function at very 
large distance. This is the behavior exhibited by 
the sequences (7.4), (7.8), (9.1), (9.2), and (9.4). 
However, there is also an approach to a limit for 
sequences like 


C(+:|r2), C(+irers| ra), 
7541), 


These probabilities show the ordering influence 
at a site immediately adjacent to stretches of 
ordered atoms of different lengths. The correla- 
tion appears to increase as the length of the 
ordered stretch increases. For stretches of six 
or more atoms, however, the correlation prob- 
abilities are identical, to terms of the order 6". 
We can conclude, then, that to this order of 
accuracy the correlation probability at a point 
immediately adjacent to a completely ordered 
half-line, r_.- is given by 


(9.5) 


In the same way, the correlation probabilities 
for sites 2, 3, 4, 5, --- steps away from the 
ordered half-line are given by 


+0|172) = — 168° 

— 648" 
-r-ito|rs) =1—6*—46*— 176° 

—756" 

t+o| 174) = — 1768 

—766" 

766" — 357g! 


(9.6) 


To terms of order 8” 


For the long range effect of the ordered half- 
line we therefore have the expression 


+0| = 46 — 176° 


— 768" —35762—---. (9.7) 


m, 
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Comparing (9.7) and (5.10) we obtain the inter- 
esting result that 


+ fea) =$(14+S). (9.8) 


The explanation of this identity would seem to 
be that the constraint of keeping an entire half- 
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line ordered provides sufficient knowledge to 
decide which sites in the lattice are a sites and 
which are 8 sites, in the sense of Section 5. The 
probability for finding a right atom would then 
be given in terms of the long range order by (9.8) 


or (5.6). 
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We have considered a two-dimensional square net consisting of four kinds of atoms supposing 
that only nearest neighbors interact and that there are only two distinct potential energies of 


interaction, one between like and one between unlike atoms. In extension of a method due to 
Onsager it is found that for the case where like atoms attract one another a simple “reciprocity” 
relation exists between the partition functions at pairs of temperatures ‘‘reciprocally” related 
to one another. As one temperature T tends to zero, the other 7* tends to infinity. If one 
further assumes that only one ‘“‘Curie”’ transition point exists, the relation between T and 7* 
enables one to locate the Curie temperature. Predictions can be made concerning the nature 
of the transition point with results similar to those of Kramers and Wannier. The reciprocity 
relation for the case of attraction between like atoms is found to be not valid for the case where 


‘unlike atoms attract one another. 


INTRODUCTION 


N a recent paper' on the statistics of two-dimen- 
sional ferromagnets, Kramers and Wannier 
have given a treatment of a model for cooperative 
phenomena. They discovered an important prop- 
erty of their system which expressed itself in the 
form of a simple symmetrical relation between 
the partition functions Z and Z* at pairs of 
temperatures T and 7%, related in a certain 
symmetrical way. Any knowledge of the partition 
function Z at temperature T implies an equiva- 
lent amount of knowledge about Z* at the related 


temperature 7*. The relation between T and 7* 


isaone toone relation; T decreases monotonically 
if 7* increases so that 7-0 as T*->~. There is 
only one temperature at which T= 7*. Assuming 


or the oO or of Phi y in the Faculty o 
Pure Science, Columbia University, Rew York, New York. 
Publication assisted by the Ernest Kempton Adams Fund 


for Physical Research of Columbia University. 
1H. A. and G. H. Wannier, Phys. Rev. 60, 
252 (1941). 


long range order at 7=0 and absence of long 
range order at T=, one will expect a Curie 
point. If the Curie point is the only temperature 
at which Z becomes singular, this temperature 
must be the one at which T and 7* become equal, 
Furthermore, one finds that the types of singu- 
larities which might occur at this transition are 
restricted by the relation between Z and Z*. For 
example, if one assumes that there is no latent 
heat at the Curie point, the specific heat must be 
either continuous through the transition or infi- 
nite on both sides of the transition temperature. 

Onsager’ has been able to obtain the reciprocity 
relation discussed above from a more direct con- 
sideration of the form of the partition function. 
Moreover, his ingenious methods appear to lend 
themselves readily to generalization. Onsager 
proceeds by giving two formally different ex- 
pressions for the partition function at a tempera- 


appeared in print. We are 


* This method has not 
sager for his permission to 
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ture J and at a temperature 7*. He is then able 
to produce a transformation between T and 7* 
which brings the first form of Z(T) into a simple 
relation with the second form of Z(7*). 

In this paper, Onsager’s methods have been 
extended to the study of two-dimensional square 
lattices containing four kinds of atoms in equal 
abundance. In addition to the assumption that 
only the interaction of nearest neighbors are im- 
portant, the further assumption is made that 
there are only two distinct energies of interaction, 
one between like and the other between unlike 
atoms. It must not be supposed that this re- 
striction reduces the problem to that of a crystal 
with only two kinds of atoms, for in the four- 
component lattice the a priori probability of 
finding a given kind of atom at a lattice point is } 
and not }. 

The detailed investigation will show that for 
the case where like atoms attract one another a 
simple symmetry or ‘‘reciprocity”’ relation, in the 
sense of Kramers and Wannier, exists. The rela- 
tion does not exist, however, if like atoms repel 
each other. One can understand this result by 
using the following qualitative consideration. At 
low temperatures, for the case where like atoms 
attract one another, minimum energy is attained 
when the atoms separate into regions consisting 
exclusively of like atoms; at high temperatures 
the atoms should be randomly distributed on the 
lattice sites. One might, therefore, expect a 
transition point where the long range order 
disappears. The situation is different when unlike 
atoms attract one another. Using the letters 
A, B, C, D to denote the four kinds of atoms, if, 
as is assumed, the unlike pairs, A~B, A—C, A—D, 
etc., attract and have the same energy of inter- 
action, there are many different states of mini- 
mum energy, and from the presence of an A atom 
on a given site, one can draw no conclusion about 
the occupant of a distant site. Therefore, there is 
no long range order at low temperatures and no 
reason to expect an order-disorder transition. 


1. FIRST FORM OF THE PARTITION FUNCTION 


Consider a two-dimensional square lattice con- 
sisting of four kinds of atoms A, B, C, D. Suppose 
that only nearest neighbors interact and that the 
energies of interaction between the various pairs 
of nearest neighbors are given by the following 


scheme: 


Pairs A-B A-C A-D 


C-D B-D B-C 


Energy € 


The reason for splitting the six pairs of unlike 
atoms into three groups of two equivalent pairs 
each will appear later. Eventually, ¢’, will 
be taken as equal. 

The first form of the partition function is given 
in terms of four auxiliary functions of the temper- 
atures u, x, y, and z defined by the relations 


exp (—¢/kT) =u+x+y+2, 
exp (—¢'/kT) =u—x+y~—z, 
exp (—e”/kT) =u—x—y+z, 
exp (—¢’"/kT) =u+x—y—z. 


By a “‘configuration”’, denoted by the general 
letter c, we shall mean any particular distribution 
of the A, B, C, and D atoms on the various 
available lattice sites. The partition function Z, 
by using the grand canonical ensemble, is defined 


by 
exp (—E./kT), 


where E, is the total energy associated with the 
configuration c, and the summation is taken over 
all possible configurations. Since E, is a certain 
linear combination of «, ¢’, with integral 
coefficients, Z can also be written as 


where the numbers of factors in the various prod- 
ucts depend, of course, on the configuration c. 
Since each pair of adjacent lattice sites con- 
tributes a factor to one of the four products in 
(2), we have, recalling the definitions of u, x, y, 
and 


(1) 


(3) 


Z=>{ (wtx+y+s)}, 
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Fic. 1. An ineffective combination of lines. 


where the product for each configuration is taken 
over all pairs of adjacent sites. This product can 
be evaluated by the following procedure: (i) for 
each pair of adjacent lattice sites, we choose a 
number which is either u, x, y, or z; (ii) examining 
the pairs of atoms assigned by the configuration 
c to the pairs of adjacent sites we attach the 
proper + or — signs to the x, y, or 2 factors 
chosen in (i) and multiply together the u, +x, 
+y, +2 factors for this choice; (iii) we repeat the 
procedure of (ii) for all possible choices made in 
(i) and add the results together, the sum giving 
the required product in (3) for the given con- 
figuration c. 

The choice indicated in (i) can be specified in a 
more geometrical way if we draw in the lattice 
connecting lines labeled with an x, y, or z between 
all pairs of adjacent lattice sites for which an x, y, 
or 2 is selected. No connecting line is drawn be- 
tween a pair for which the letter u is chosen in (i). 
In this way the choice in (i) is equivalent to a 
“pattern” of x, y, and z lines running through the 
lattice. Any configuration c of atoms in the 
lattice may be adopted for any pattern of lines p. 
For a given pattern p, the choice of signs in (ii) is 
determined by the configuration c. Therefore, the 
partition function may be written in the form of 
a double sum, summations being carried out over 
all configurations c, and over all patterns p. 

ce patterns 
The number of factors in each of the four products 
depends on the pattern p alone. The only effect 
of the configuration c is to influence the distri- 
bution of + and — signs. 

There are many patterns, however, which do 
not effectively contribute to Z. In investigating 
the combinations of x, y, and z lines of a pattern 
which can intersect at a lattice point, we find 
that, under certain circumstances, there is a 
mutual cancellation of the four contributions to 
Z arising from a given pattern and four configu- 
rations which differ among themselves only in the 
identity of the atom at the lattice intersection 
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point in question. For example, suppose that at a 
point a single x line and a single y line terminate. 
For any configuration, the contribution to Z of 
a pattern containing this combination of x and y 
lines intersecting, will evidently have a factor 
(+x) (+y) arising from these intersecting lines, 
Let us consider this factor for each of the four 
configurations obtained from some arbitrary 
initial one by successively changing the identity 
of the atom V at the point of intersection, keeping 
the identities of all other atoms fixed. Looking at 
the definition of u, x, y, z for the signs, we have 
(Fig. 1) an ineffective combination of lines. The 
sum of these contributions vanishes. The same 
result is found for other initial configurations 
which place other atoms at the points 1 and 2, 
We conclude, therefore, that any pattern which 
contains the combination x y at a point, does not 
effectively contribute to Z. 

A different result (Fig. 2) is obtained for the 
combination of intersecting lines x y z. All these 
factors are of the same sign; their total contribu. 
tion does not vanish. The same conclusion is 
reached for other atoms at points 1, 2, and 3. 

Examining in this manner all the various com- 
binations of lines which can intersect at a point, 
we find that for those patterns which are effective 
for Z, ‘‘effective patterns” for short, the combi- 
nations of intersecting lines must be among the 
following types: 


xx xyz xxXXX xxyy 


vo 
23 2222 yy ze. 


Figure 3 gives an example of a typical effective 
pattern. 

The effective combinations listed above have 
the property that the signs of their contributions 
to Z remain unchanged however one changes the 
identity of the atom at the point of intersection. 
From this we can prove that the contributions to 
Z of a given effective pattern p, are positive for 
all configurations of atoms ¢, and therefore are all 
equal and independent of c. This is evident if we 


ATOMV| A 8 C 


Fic. 2. An effective combination of lines. 
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consider that the contribution for a configuration 
of A atoms on all lattice points is surely positive, 
and that successive changes in the identities of 
the atoms at separate lattice points preserve the 
sign. 

If N is the number of lattice sites, the total 
number of configurations is 4%. Confining the 
summation to effective patterns and using the 
result of the previous paragraph, we can write as 
the final expression for the first form of the 
partition function 


eff. patt. 


(5) 


where /, m, n, are the total lengths (in units of the 
lattice distance) of the x, y, z lines, respectively, 
in one of the effective patterns contributing to 
the sum. Since there are 2N pairs of nearest 
neighbors when WN is a large number and the 
factor u is chosen for all pairs of adjacent sites 
not included in the pattern, u appears to the 
power 2N—(/+m-+n). 


2. SECOND FORM OF THE PARTITION 
FUNCTION 


The second form of the partition function is 
obtained by a geometrical procedure different 
from the one just used. For a given configuration, 
let us pick out all pairs of unlike nearest neigh- 
bors. Then, for each such pair, let us draw a 
separating line midway between the two atoms, 
this separating line being perpendicular to the 
line joining the pair and of length equal to the 
lattice distance. These separating lines will, from 
their definition, form the boundaries between 
regions in the lattice consisting of A atoms alone, 
B atoms alone, etc. 


Fic. 3. A typical effective pattern. 


Fic. 4 (left). An impossible combination of lines. 
Fic. 5 (right). A possible combination of lines. 


Let us define four functions of the temperature 
a, B, y, 6 by 


€ € 


and proceed to label the boundary lines ac- 
cording to the following scheme in accordance 
with the fact that the energy between an A—B or 
C-D pair is ¢’ and so on. 


Pairs of A-B A-C A-D 
unlike atoms C-D B-D B-C 

Label of B Y é. 
separating line 


In this way, each configuration has associated 
with it a pattern of 8, y, 6 lines. These patterns 
are not entirely arbitrary, however. For example, 
we cannot have a single 8 line and a single y line 
terminating at a point. In Fig. 4, a B atom was 
assumed at 1, and the atoms 2 and 3 written 
down in accordance with our labeling scheme. 
The result is a contradiction, since the atoms at 
2 and 3 are not separated by any line and must 
therefore have the same identity.* On the other 
hand, the combination 8 y 4 is possible, as Fig. 5 
shows. If the identity of atom 1 is specified, the 
identities of atoms 2, 3, and 4 follow uniquely. 

Examining all the cases, we find that the only 
possible combinations of 8, y, 6 lines which can 
intersect at a point are among the following 


types: 


BB BBBB 
YY YYYY 
66 6666 


BByy 


* The argument holds only for the interior of the two- 
dimensional crystal. The resulting reciprocity relation is 
proved, therefore, only for the case where surface effects 
may be neglected. 
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L_ 


é 


4 8 8B 

Cc 


Fic. 6. 8, y, 6 pattern identical with the x, y, z 
pattern of Fig. 3. 


Cc 


From this it can be concluded that, except for 
effects at the edges of the lattice which can be 
neglected when the lattice is very large, the 
patterns of 8, y, 6 lines are geometrically identical 
with the effective patterns of x, y, z lines of the 
previous section. Figure 6 shows a pattern of 
B, y, 6 lines identical with the pattern of x, y, z 
lines of Fig. 3. 

It is evident that each effective, or allowed, 
_ pattern of 8, y, 6 lines is associated with exactly 
four distinct configurations. For if we give the 
identity of a particular atom in the lattice, the 
identities of all other atoms follow uniquely from 
the 8, y, 6 pattern. The four configurations arise 
from the four possible identities of this particular 
atom. The energy of each of these configurations 
is 


(7) 


where /, m, n are the respective lengths (in units 
of the lattice distance) of the 8, y, 6 lines of the 
associated pattern. Remembering the factor four 
and the definitions of a, 8, y, 6, we can write the 
partition function as a sum over patterns, instead 
of configurations, giving 
Z=4 
eff. patt. 
The next step is to compare the expressions (5) 
and (8). 


3. THE RECIPROCITY RELATION 


We now introduce the assumption that all 
pairs of unlike nearest neighbors interact with 


(8) 


the same energy; that is, we take e’=¢’’=e’’’. As 
a consequence x=y=z and B=y=6. For con- 
venience we shall write the energies ¢ and ¢’ in © 
the form 


e=—wt¢ 
and instead of the temperature, introduce as the 


new independent variable the quantity 6 defined 
by 


(9) 


6=exp (w/kT)=aexp (¢/kT). (10a) 


We also define three new functions of the temper- 
ature U, X, 7 


U=uexp ($/kT), 
X =x exp (¢/kT), 
n=1/0=8 exp (¢/kT). 


Then from (1), (9), (10) and the equality of x, y, 
and 2, it follows that 


exp [(w—¢)/kT ]=u+3x 
=exp (—¢/kT)(U+3X), 


exp [(—w—$)/kT]=u—x 
=exp (—$/kT)(U—X). 


The quantity ¢ can evidently be eliminated from 
these equations and one finds for U and X as 
functions of 6, the expressions 


X= 7[6—(1/6)], 
U=}[0+ (3/6) 


Let us now consider the partition function at 
two temperatures T and 7* which are, so far, 
unrelated to one another, writing Z = Z(T) in the 
form (5) and Z*=Z(7*) in the form (8). For 
1+m-+n we shall write ¢. Using (10) to express 
u, x, a, B in terms of U, X, 6, 7 we have 


Z=4% u2N—txt 
eff. patt. 


=4% exp(—2N¢/kT)U* (X/U)', (13) 


eff. patt. 
Zt=4 a*2N—tget 
eff. patt. 


=4exp (n*/6*)! 


eff. patt. 


=4exp(—2No/kT*)0*2% (1/0**)!. (14) 


eff. patt. 


(10b) 


(12) 
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STATISTICS OF TWO-DIMENSIONAL LATTICES 


By a “reciprocity” relation between Z and Z* 
we understand the following. There shall exist a 
transformation 7*= F(T) between the tempera- 
tures J and 7* having the two properties (i) if T 
is carried into J7* by the transformation, the 
transformation also carries T* into 7; that is, if 
T*=F(T), then T= F(T*), (ii) for this transfor- 
mation there shall exist a simple function G(T) 
whose form is readily obtained without explicity 
evaluating the partition function and which is 
such that 


(18) 


Moreover, with the possible exception of the 
points T=0 and T=~, G(T) shall have no 
singularities for real temperatures. If these con- 
ditions are satisfied one can conclude that if 
Z(T) has a singularity at T, Z(T*) will have a 
corresponding singularity at T*. 

A comparison between (13) and (14) suggests 
that the transformation between T and 7* is the 
one given implicitly in the condition 


= U/X. (16) 


We shall verify that (16) leads to a reciprocity 
relation in the sense previously described. 

Writing U and X in terms of 6 we see that (16) 
is equivalent to 


= (6?+3)/(6?—1). 
Solving for @ gives 
6? = (18) 


so that the transformation from T into 7* is of 
the required type. 

Substituting (16) into (14) and comparing with 
the expression for Z in (13) we have 
{exp (2Ng/kT*)Z* = (U/X)* (X/U)! 


=(4UX)- exp (2N@/kT)Z. 


(17) 


(19) 


The relation to be derived is most simply ex- 
pressed in terms of the Nth root of Z, hereafter 
denoted by \. This quantity may be regarded as 
the partition function per particle. In taking this 
Nth root, however, we shall replace 4" by unity. 
To this extent, and to the extent that effects due 
to the edge of the crystal are neglected, our 
method of deriving the reciprocity relation is 
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restricted to very large crystals. Writing \=Z"/* 
and performing the Nth root we have 


exp (2¢/kT*)A* = (4UX)— exp (2¢/kT)A. (20) 


Using (12), (17), and (20) one may verify the 
expression 


exp (26/kT)A = exp (26/kT*)A* (21) 


obtained in a formal way from (20) by inter- 
changing starred and unstarred quantities. From 
(20) and (21) the symmetrical relation 


exp (2¢/kT*)\* 
= exp (2¢/kT)A 


readily follows. This has the form of a reciprocity 
relation. 
If we define functions f(@) and ¥(@) by 


f(0) 
¥(6) =f(8)A(8) 
the reciprocity is also expressed by 
¥(6*) =y(6). (24) 


The f(6) defined here plays the role of G(T) 
of (15). Since its only singularities are at @= © 
and @=1 or and T= ~, this requirement on 
G(T) is also satisfied. 

In the formal procedure outlined so far it has 
not been specified whether like or unlike atoms 
attract one another. For the case where like 
atoms attract one another, we have ¢’ greater 
than e and hence w>0. The quantity @ is therefore 
greater than unity and from (17) 6* is positive 
and also greater than unity. Thus the transforma- 
tion (17) is entirely consistent and the reciprocity 
relation (22) is valid. 

- However, for the case where like atoms repel 
one another, the energy e’ must be less than the 
energy €, or w must be negative. Consequently 
6<1. From (17), 6** would then be negative and 
would not correspond to any real temperature 7™*. 
The conclusion is, therefore, that for repulsion 
between like atoms no reciprocity relation, in the 
sense defined above, is found. 

It has been pointed out in the introduction that 
long range order at low temperatures exists if like 
atoms attract but not when they repel each other. 
Therefore, in these two cases, an essentially 


(22) 
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different behavior of the partition function at 
low temperatures is to be expected while no such 
difference at high temperatures is anticipated. 
Thus it is reasonable to find that a correlation 
between behaviors at high and low temperatures 
found for the one case does not hold for the other. 
It must not be supposed that no relation of the 
form (15) can exist for the case of repulsion be- 
tween like atoms. In fact an infinite number of 
such relations is possible. As an example, suppose 
T and 7* are related by TT7T*=T;,. Let 
G(T)=Z(T¢/T). Then we can evidently write 


G(T)Z(T) =G(T*)Z(T*). 


But now no conclusion about Z(7) can be drawn 
since the form of G(T) remains unknown as long 
as that of Z(T) isnot known. 

Returning to the case of attraction between 
like atoms we see that from (17) 7* is a mono- 
tonically decreasing function of T. As T increases 
from 0 to «, 7* decreases from ~ to 0. There is 
only one temperature, which we shall denote by 
T., for which T= 7*. From (17) this occurs at 


6.=exp (w/kT.) = V3. (25) 


The interval from 0 to 7, is mapped by the 
transformation upon the interval from to 
Knowledge of \ on one of these intervals is 
equivalent, by (22), to knowledge of it on the 
other. If \ has any singularity at temperature 7, 
it will have the same type of singularity at T*. If 
it is assumed that only one singularity exists, it 
must then occur at temperature 7., which we 
shall now call the critical temperature. 

In Eq. (25) the quantity w is half the difference 
between the energies of interaction between like 
and unlike atoms. The critical point for the two 
component system of Kramers and Wannier is 


given by 
0.2 =exp (2w/kT.) =1+Vv2, (26) 


where w has the same significance as in our case. 
It will be noted that, for the same numerical 
value of w, the transition occurs in our case at a 
temperature lower than the transition tempera. 
ture of Kramers and Wannier by the ratio of 
In (1+2) to In3 or by a factor of about 0.8, 
This may be due to the fact that in the four. 
component system the a priori probability favors 
disorder more than in the two-component system 
so that long range order breaks down at a lower 
temperature. This consideration gives support to 
the view that @, corresponds to the ‘‘Curie point.” 

The existence of the reciprocity relation (22) or 
(24) enables us to make certain limited state. 
ments concerning the nature of the hypothetical 
transition point. These conclusions are identical 
with those of Kramers and Wannier for lattices 
with two components. The energy E and specific 
heat C are given in terms of \ by 


E=-—Nd log \/d(1/kT) = — Nw6d log d/d8, 
C=dE/dT =(—w0/kT*)(dE/d8). 


From (17) and (24) we find for the possible 
discontinuities in E and C at the critical point 


AE= —2Nw0.(d log +0, 
AC= (w/kT.*) AE. 


Thus, if the energy is discontinuous, so is the 
specific heat. If AE=0 the specific heat is either 
continuous or becomes infinite at the transition 
temperature in such a way that the difference in 
its values at T and 7* tends to zero as T tends 
to T.. This is about as far as one can go with the 
reciprocity relation alone. 


(27) 
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ROMPT publication of brief reports of important dis- 

coveries in physics may be secured by addressing them 
to this department. The closing date for this department is the 
third of the month. Because of the late closing date for the sec- 
tion no proof can be shown to authors. The Board of Editors 
does not hold itself responsible for the opinions expressed by 
the correspondents. Communications should not in general ex- 
ceed 600 words in length. 


The Shape of Betatron Pole Faces 


James H. BARTLETT 
Department of Physics, University of Illinois, Urbana, Illinois 
August 25, 1943 
OR there to be focusing of the electron beam in the 
betatron, it is necessary that the magnetic field near 
the beam fulfill certain requirements.’ In practice, these 
have been met by selecting a certain shape of pole face, 
determining the field by approximate methods such as 
flux plotting, and then modifying the original shape if it 
was not a suitable one. It is the purpose of this note to 
point out that the pole face shape may be determined in a 
more direct manner. 

As in the article of Kerst and Serber, we use cylindrical 
coordinates, and z= 0 describes the median plane. In empty 
space, curl H=0 and div H=0, so that H=—grad ¢n, 
where ym is the magnetostatic potential, and y*¢,, =0. We 
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shall assume that H,(r) at z=0 is a known function, which 
is such that it satisfies the focusing conditions. A surface 
of constant magnetostatic potential will be, to a high degree 
of approximation, an acceptable boundary for the magnet. 
The problem thus reduces to the calculation of ¢, subject 
to the condition at z=0. In what follows, we drop the 
subscript m. 
We require the solution of: 
1 


or 
Expand ¢ as a power series in z, so that: 


o= (1) 
n=0 


Substituting in the differential equation, one obtains the 
recursion formula: 


Now: 

H,= — (8¢/dz) = — Zn 
Therefore, (H,)z.0= — ¢1. The potential go is arbitrary and 
may therefore be set equal to zero. The solution (1) of the 
problem is thus obtained by differentiating expressions in- 
volving ¢: [as prescribed by (2) ], and will be an odd func- 
tion of s. 

If the series (1) should prove to be valid in too narrow a 
range, one could extend the range by analytic continuation, 
i.e., a power series in —20, where Zo is such a value that 
¢ may be readily calculated by (1). 


1 Stated by Kerst and Serber, Phys. Rev. 60, 58 (1941). 
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PHYSICAL REVIEW 


HE 256th meeting of the American Physical 

Society was held at Stanford University, 
on Saturday, July 10, 1943. In spite of the 
“unusual” conditions prevailing on the campus, 
due to the presence of very large numbers of 
students in the armed forces, completely satis- 
factory arrangements for the meeting had been 
made, in accordance with the usual Stanford 
tradition. The average attendance at the two 
sessions was about 55, which nearly equals the 
average attendance at previous Stanford meet- 
ings. Some fifty members of the society and 
guests had luncheon together in the banquet 
room of Wilson’s Restaurant, Palo Alto. Brief 
remarks were made at the luncheon by Professor 
Paul Kirkpatrick who, as chairman of the host 
department, also presided at the scientific ses- 
sions; by Dr. K. K. Darrow; and by the Local 
Secretary. The chief topic of these remarks was 
the trying situation in which all physics depart- 
ments now find themselves, as the result of the 
specialized programs now being offered to large 
numbers of the armed forces. 


1. The Validity of Lens Equation and Magnification 
Formula of Light Optics for Electron Lenses. R. G. E. 
Hutter, Division of Electron Optics, Stanford University.— 
The usefulness of the lens equation (Zo—ZFo) -(Z;—ZF;) 
= — F; of the magnification formula M = Fo/(Zo—ZFo) 
of light optics, and of the simple geometrical image con- 
struction based upon these equations, lies in the fact that 
the set of lens constants-focal points ZFo, ZF;, and focal 
lengths Fo, F; are characteristics of the lens and are inde- 
pendent of the object position. It has often been shown that 
these equations also apply to electron optical lenses. These 
proofs assume, however, that the lens fields are of finite 
extent on the optical axis and that object and image are 
outside the lens region. The assumptions are of the nature 
of sufficient conditions only; i.e., electron lenses exist for 
which none of these conditions holds and yet the lenses 
can be characterized by a single set of optical constants. 
Necessary and sufficient conditions for such lens types are 
stated and an analytical expression for all lenses is derived. 
A number of simple magnetic and electrostatic lenses 
satisfying these conditions are discussed. 


VOLUME 64, NUMBERS 5 AND 6 
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Because of an expected dearth of contributed 
papers, arrangements had been made for three 
invited papers, and the success of the meeting 
was undoubtedly due in great measure to these 
papers. They were as follows: 


Recent Developments in Electron Microscopy. L. Mar.- 
TON, Division of Electron Optics, Stanford University. , 

The Applications of Nuclear Physics to Medicine and 
Biology. HARDIN JONES, Radiation Laboratory, University 
of California. 

Entropy and Probability. K. K. Darrow, Bell Telephone 
Laboratories. 

Dr. Hardin Jones, in connection with his paper, 
gave experimental demonstrations with an auto- 
matic recording Geiger counter, using a live 
mouse, on which he performed a surgical opera- 
tion in full view of the audience. 

There were twelve contributed papers, for 
which the abstracts are appended. None of these 
was read by title, which, I believe, constitutes a 
new record. 

R. T. BirGE 
Local Secretary for the Pacific Coast 


2. Optimum Conditions for Apertures of Magnetic Elec- 
tron Lenses of the Field Form H(Z) = Ho/(1+(Z/a)®). L. 
MARTON AND R, G. E. Hutter, Division of Electron Optics, 
Stanford University——It has been shown previously that 
the optical properties of strong magnetic lenses of the 
field form H(Z) = Ho/(1+(Z/a)*) can be completely calcu- 
lated. Such calculations up to now did not include the 


position of the lens aperture. Furthermore, the size of the’ 


aperture was derived from considerations of the aberrations 
of weak lenses only. By using the complete analytical 
expression for the electron path through a lens, both loca- 
tion and size of such an aperture can be calculated. This 
aperture achieves the ray-limiting effect at a maximum 
value of its diameter. The calculated size and location of 
the aperture differs considerably from the values currently 
adopted. Diagrams are given for the aperture size and 
location in function of such parameters which allow 
their application to lenses of different strength and mag- 
nification. 
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3. Intensity of High Energy Electrons and Photons at 
10,000 Feet. Wayne E. Hazen, University of California.— 
A cloud chamber containing eight 0.7-cm lead plates was 
operated without counter control at an altitude of 10,000 
ft. High energy electrons and photons produced cascade 
showers in the lead plates whereas mesotrons produced 
only an occasional collision electron. Usually the shower 
maxima were observed, since the 30 by 30 cm chamber size 
provided a large solid angle and there was sufficient lead 
to produce maxima even in showers of several hundreds of 
particles. In 8500 photographs, 1090 showers that reached 
a maximum of four or more particles were observed. If 
we use the number of particles at the maximum as a 
measure of the energy of the incident particle, the differ- 
ential energy spectrum for electrons plus photons of 
energy 2X 10* to 10° ev follows an E~*+ distribution. When 
combined with Greisen’s' counter data for electrons with 
energy greater than 2X10’ ev, the present results indicate 
an E-*-* distribution. The intensity of electrons plus 
photons with energies greater than 2X 10* ev relative to 
penetrating particles was found to be 8.5 percent. 


1K. Greisen, Phys. Rev. 61, 212 (1942). 


4. A Possible New Application of the Wilson Cloud 
Chamber. Cart E. Nietsen, University of California.— 
It is well known that the atmosphere contains objects that 
serve as condensation nuclei with very small supersatura- 
tions. Such nuclei—which have interested physicists 
primarily in being objects that must be removed before 
the Wilson cloud chamber is useful for the study of ions— 
may, it is commonly supposed, be inert dust particles as 
well as salt particles and certain other chemicals. However, 
it has been asserted that condensation cannot occur upon 
“ordinary” dust particles. The experiments to be described 
have shown that dust (carbon) particles 0.5 micron in 
diameter, and also bacteria of dimensions 0.5 micron 
to 2 microns do serve as condensation nuclei with 
small supersaturations. Consequently it is possible to use 
the cloud chamber to precipitate bacteria from the air 
upon a petri dish containing agar culture medium. Unlike 
most common air sampling techniques, which give an 
unknown percent removal, this method removes all of the 
bacteria from the air. 


5. Townsend Currents in Non-Uniform Electric Fields. 
Paut L. Morton, University of California. (Introduced by 
Leonard B. Loeb.)—The current in the Townsend discharge 
in nearly uniform fields and in the absence of secondary 
effects can be calculated from the equation log i/ip= fadx, 
where the coefficient a as a function of the field intensity 
can be taken from the experimental work of Hale and 
others. Uniform fields seldom occur because of space 
charge distortion as well as electrode shape. The accuracy 
of this calculation has therefore been investigated for the 
field between concentric cylinders, in hydrogen at pres- 
sures from 0.01 to 10 millimeters, with currents sufficiently 
small to avoid space charge effects. The results indicate 
that large errors should be expected when the field inten- 
sity near the cathode changes by more than 10 percent 
per electron mean free path even when the gap is com- 


paratively long, and that the calculated currents can be 
orders of magnitude off when the electrons reach the 
anode before attaining a terminal drift velocity corre- 
sponding to the field there. If the field changes rapidly 
only near the cathode and the gap is long, the total ioniza- 
tion as calculated may be approximately correct, even 
when its location, as indicated by the magnitude of a, is 
greatly in error. This fact accounts for the approximate 
agreement sometimes found when the Townsend equation 


is applied to glow discharges. 


6. Apparent Failure of the Meek Criterion for Streamer 
Formation. Leon H. FisHer, University of California.— 
The criterion for streamer formation given by Meek is 
that the positive space charge field of an electron avalanche 
be equal to K times the impressed field. Measurements of 
streamer onset in confocal paraboloid gaps by Weissler 
and Loeb! indicate that K is extremely pressure dependent. 
In hydrogen and dry air, K decreases rapidly with de- 
creasing pressure. The possibility existed that the method 
used in calculating the space charge field in a non-uniform 
gap was not valid. To test this, the cathode of the least 
uniform gap of Weissler and Loeb was illuminated with 
ultraviolet light, and measurements of the current below 
onset were made in hydrogen. The results indicate that 
the ionization in these divergent fields cannot be calculated 
by integrating the ionization functions for uniform fields. 
At a constant pressure, the current-voltage curve increases 
too rapidly with voltage. This effect becomes more pro- 
nounced with decreasing pressure. Morton has obtained 
similar results (paper No. 5) in much less divergent fields 
in hydrogen at pressures near one millimeter. The present 
results explain the observed trend in K, and restore the 
possibility of its physical significance. Measurements in air 
were also attempted, but no currents below streamer onset 
could be measured. This observation can be explained by 
electron attachment to oxygen molecules in the weak 
portion of the field. ; 


1 Weissler and Loeb, Phys. Rev. 62, 300 (1942). 


7. The Effect of Adsorbed Gas on the Contact Electrifi- 
cation of Quartz on Nickel. Davip E. DeBeau, University 
of California.—The contact electrification of quartz par- 
ticles in falling from one Ni container to a second Ni con- 
tainer has been measured as a function of gas pressure. 
When the gas used is air or O2, the magnitude of the charge 
collected is a minimum at pressures of about 0.1 to 1.0 mm 
of Hg. The charge collected at 10-§ mm and at 760 mm is 
approximately the same. The maximum charge collected 
in one run is about 80 e.s.u. and the quartz and the Ni 
receiver become negative. When H; is used, no consistent 
results are obtainable and the charge collected is low. 
Not until the chamber is baked out in the presence of air 
are consistent results again obtainable with air or O; in the 
chamber. 


8. Time Variation of Fog and Pressure in a Wilson Cloud 
Chamber. Morris Neustapr, University of California.— 
Fog and pressure were measured electrically, and were 
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observed as traces on an oscilloscope, an electron switch 
showing them simultaneously. The chamber was operated 
(1) saturated with water, (2) dry, (3) saturated with alco- 
hol-water solution, 3 to 1 by volume. All measurements 
were made both for expansion to constant final volume, 
and. to constant final pressure. Expansion ratios ranged 
from below the ion limit to above the fog limit. The larger 
the expansion ratio the larger is dp/dt after expansion. 
When large fog is present there is a rapid build-up of 
pressure due to heat of condensation of the fog. This effect 
on the pressure is not present when a dry chamber is used. 
If one expands below the ion limit, dense fog present from a 
previous over-expansion causes the pressure minimum to 
be higher, and causes the pressure to rise very sharply for 
a short time. As R increases, the fog builds up sooner and 
has a larger maximum value. Initially the fog increases as 
#, but during the sensitive time this relationship gradually 
changes, and the fog increases as ¢ after the end of the 
sensitive time. Constant pressure fogs scatter more light 
than constant volume fogs. ' 


9. Spectrographic Determination of Lead in Pectinous 
Materials. C. RULON JEPPESEN* AND E. JoHN EASTMOND, 
Western Regional Research Laboratory.—A spectrographic 
method of analysis for lead in pectinous materials has 
been developed which utilizes dry-ashed samples and thus 
avoids possible contamination from chemical treatments 
and the difficulties encountered in putting certain of these 
materials into solution. The dry-ashed sample is mixed with 
four times its weight of a standard buffer base consisting 
of 0.175 percent bismuth in lithium carbonate. A 250-volt 
d.c. arc between graphite electrodes is used as a source of 
excitation. The ratio of the intensities of the lines Pb 2833A 
and Bi 2898A is determined from photometric measure- 
ments of the spectrogram and the concentration of lead 
is read from a previously prepared working curve. Con- 
centrations of lead as low as 0.01 part per million in the 
original sample and 60 parts per million in the ash have 
been measured. Repeat analyses on the same sample show 
an average deviation from the mean of approximately 15 
percent when triplicate exposures are used. Tests indicate 
variations are due chiefly to the ashing and mixing proc- 
esses. The method is rapid and entirely objective. 

* Now at Montana State University, Missoula, Montana. 


10. Contribution to the Theory of Liquids. AHLBORN 
WHEELER, Shell Development Company, Emeryville, Cali- 
fornia. (Introduced by O. Beeck.)—To date all theories of 
the liquid state have been seriously handicapped by the 
fact that we cannot assign the correct statistical weight 
(entropy) to a given arrangement (energy state) of the 
molecules in our system. It is shown how this can be done 
for the simplest type of liquid system: a linear chain of 
molecules with arbitrary interaction potential between 
nearest neighbors. The exact partition function is con- 
structed for this system and the thermodynamic properties 
are calculated for potential functions of the Lennard-Jones 
type. The nucleus of the solution is an exponential integral 
of (E+pV)/kT over the coordinate space, and this integral 
gives the Gibbs free energy directly. The results should be 
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of considerable assistance in formulating improved theories 
for real liquids. First of all the form of the solution is very 
suggestive in indicating how to construct new and perhaps 
more accurate partition functions for the three-dimensional 
case. Secondly we have a ‘‘yardstick” with which to test the 
validity of any approximate model of a liquid. This mode} 
may be applied to our one-dimensional system, and the 
results compared with the exact solution. 


11. Identification of Hydrocarbons by Means of Their 
Infra-Red Spectra. R. S. Rasmussen, R. R. Bratram, 
AND O. BEECK, Shell Development Company, Emeryville, 
California.—Infra-red absorption spectroscopy in the 
rocksalt region (2-154) has been applied to the rapid 
qualitative and semi-quantitative analysis of higher iso. 
meric paraffin and olefin hydrocarbons (C¢ to Cs). Spectra 
of pure hydrocarbons serve as comparisons, against which 
the spectra of new preparations are checked. Samples which 
are nearly pure are readily identified in this way, and the 
type and amount of impurity can be estimated. Since ip © 
complex mixtures the bands of the individual hydrocarbons 
often overlap to a great extent, a rough distillation may be 
necessary in order to divide the mixture into various cuts 
which in turn can be examined by the infra-red absorption 
method. Further use of infra-red spectra can be made in 
those cases where characteristic bands exist for specific 
groups. These bands may be used in establishing the struc- 
ture of molecules. While this procedure is uncertain in the 
case of paraffin hydrocarbons, it has proved highly useful 
in the study of olefins. Thus the C= CH group (a-position 
of the double bond) is indicated by a band at 3.2y anda 
strong band in the 6u region. Moreover, several very 
strong bands in the 10-12.5u region are highly character- 
istic of the type of substitution around the double bond. 


12. Quantitative Analysis of Multicomponent Miz- 
tures by Means of Infra-Red Spectrophotometry. R. R. 
Brattain, R. S. RAsMuSSEN, AND A. M. Cravatu, Shell 
Development Company, Emeryville, California.—Provided 
sufficient differences exist in the absorption spectra of the 
pure components and provided Beer’s law is obeyed, spec- 
trophotometric analysis of a multicomponent mixture of # 
components can be obtained by making optical density 
measurements at m wave-length positions, and solving s 
linear simultaneous equations in » unknowns. For routine 
work in the infra-red region from 2 to 15y, desirable stabil- 
ity is gained by using wide spectral slits and no further 
amplification of the thermocouple output. Under these 
conditions Beer’s law cannot be applied to the absorption 
measurements because of the rapid change of the true 
absorption coefficient over the spectral slit width. In 
order to analyze mixtures under these conditions, it is 
assumed that the optical density of the mixture is equal 
to the sum of the optical densities of the pure components, 
although the optical density of a component is no longer 
linear with concentration. Calibration curves can be ob 
tained from data on pure compounds giving optical density 
as a function of concentration (non-linear) and can be used 
in the analysis of mixtures. The method has been thor 
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oughly tested on a large number of hydrocarbon mixtures 
and was found to give excellent results with mixtures of 
6 and 7 components. An accuracy of +0.5 percent (of total 
sample) is obtained for each component. By using a 
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specially devised rapid method of calculation the time 
required to analyze a six-component C, hydrocarbon mix- 
ture is 1} hours as against 12 to 15 hours by low tempera- 
ture distillation and subsequent gas analyses. 
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